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Abstract. Let S be the boundary of a convex polytope of dimension d+1, or 
more generally let 5* be a convex polyhedral pseudomanifold. We prove that S has a 
polyhedral nonoverlapping unfolding into M.'^, so the metric space S is obtained from 
a closed (usually nonconvex) polyhedral ball in by identifying pairs of boundary 
faces isometrically. Our existence proof exploits geodesic flow away from a source 
point V G which is the exponential map to S from the tangent space at v. We 
characterize the cut locus (the closure of the set of points in S with more than 
one shortest path to v) as a polyhedral complex in terms of Voronoi diagrams on 
facets. Analyzing infinitesimal expansion of the wavefront consisting of points at 
constant distance from v on S produces an algorithmic method for constructing 
Voronoi diagrams in each facet, and hence the unfolding of S. The algorithm, 
for which we provide pseudocode, solves the discrete geodesic problem. Its main 
construction generalizes the source unfolding for boundaries of 3-polytopes into K.^ . 
We present conjectures concerning the number of shortest paths on the boundaries 
of convex polyhedra, and concerning continuous unfolding of convex polyhedra. We 
also comment on the intrinsic non-polynomial complexity of nonconvex manifolds. 



Date: 4 December 2003. 

Key words and phrases. Convex polyhedra, nonoverlapping unfolding, polyhedral metric, cut 
locus, discrete geodesic problem, geodesic Voronoi diagram, convex polyhedral pseudomanifold. 

* School of Mathematics, University of Minnesota, Minneapolis, MN; ezra@math.uinn.edu. 

'''Department of Mathematics, MIT, Cambridge, MA; pak@inath.mit.edu. 



Contents 



Introduction 

Overview 

Methods 

1. Geodesies in polyhedral boundaries 

2. Cut loci 

3. Polyhedral nonoverlapping unfolding 

4. The source poset 

5. Constructing source images 

6. Algorithm for source unfolding 

7. Convex polyhedral pseudomanifolds 

8. Limitations, generalizations, and history 

9. Open problems and complexity issues 
Acknowledgment s 

References 



121 
121 
il 
□ 

cni 
nil 

1221 
I2S1 
M 
EH 
iU 
M 
M 



1 



2 



EZRA MILLER AND IGOR PAK 



Introduction 

The past several decades have seen intense development in the combinatorics and ge- 
ometry of convex polytopes |Zie95j . Besides their intrinsic interest, the advances have 
been driven by apphcations to areas ranging as widely as combinatorial optimization, 
commutative algebra, symplectic geometry, theoretical physics, representation theory, 
statistics, and enumerative combinatorics. As a result, there is currently available a 
wealth of insight into (for example) algebraic invariants of the face posets of poly- 
topes; arithmetic information connected to sets of lattice points inside polytopes; 
and geometric constructions associated with linear functionals, such as Morse-like 
decompositions and methods for locating extrema. 

On the topological side, there are metric theories for polyhedral spaces, primarily 
motivated by differential geometry. In addition, there is a vast literature on general 
convexity. Nonetheless, there seems to be lacking a study of the interaction between 
the combinatorics of the boundaries of convex polytopes and their metric geometry in 
arbitrary dimension. This remains the case despite relations to a number of classical 
algorithmic problems in discrete and computational geometry. 

The realization here is that convexity and polyhedrality together impose rich com- 
binatorial structures on the collection of shortest paths in a metrized sphere. We 
initiate a systematic investigation of this metric combinatorics of convex polyhedra 
by proving the existence of polyhedral nonoverlapping unfoldings and analyzing the 
structure of the cut locus. The algorithmic aspect, which we include together with 
its complexity analysis, was for us a motivating feature of these results. That being 
said, we also show that our general methods are robust enough so that — with a few 
minor modifications — they extend to the abstract spaces we call 'convex polyhedral 
pseudomanifolds', whose sectional curvatures along low-dimensional faces are all pos- 
itive. To conclude, we propose some directions for future research, including a series 
of precise conjectures on the number of combinatorial types of shortest paths, and on 
the geometry of unfolding boundaries of polyhedra. 

Overview 

Broadly speaking, the metric geometry of boundaries of 3-dimensional polytopes is 
quite well understood, due in large part to work of A. D. Aleksandrov |Ale48t lAleSOj 
and his school. For higher dimensions, however, less theory appears in the literature, 
partly because Aleksandrov's strongest methods do not extend to higher dimension. 
Although there do exist general frameworks for dealing with metric geometry in spaces 
general enough to include boundaries of convex polyhedra, such as |BGP92j . the 
special nature of polyhedral spaces usually plays no role. 

The existing theory that does appear for polyhedral spaces is motivated from the 
perspective of Riemannian geometry, via metric geometry on simplicial complexes, 
and seems mainly due to D.A.Stone; see |Sto76j . for example. In contrast, our 
original motivation comes from two classical problems in discrete and computational 
geometry: the 'discrete geodesic problem' |MitOOj of finding shortest paths between 
points on polyhedral surfaces, and the problem of constructing nonoverlapping un- 
foldings of convex polytopes lO'Rnnj . Both probl ems are well understood for the 
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2-dimensional boundaries of 3-polytopes, but have not been attempted in higher di- 
mensions. We resolve them here in arbitrary dimension by a unified construction 
generahzing the 'source unfolding' of 3-dimensional convex polyhedra |VP7H l5S86j . 

Previous methods for source unfoldings have been specific to low dimension, relying 
for example on the fact that arcs of circles in the plane intersect polygons in finite 
sets of points. We instead use techniques based on differential geometry to obtain 
general results concerning cut loci on boundaries of polytopes in arbitrary dimen- 
sion, namely Theorem 12 . 91 and Corollarv 12.111 thereby producing polyhedral foldouts 
in Theorem 13.51 In more precise terms, our two main goals in this paper are to: 

1. describe how the set of points on the boundary S* of a convex polyhedron at given 
radius from a fixed source point changes as the radius increases continuously; 

2. use this description of 'wavefront expansion' to construct a polyhedral nonover- 
lapping unfolding of the (i-dimensional polyhedral complex S into M°'. 

By 'describe' and 'construct' we mean to achieve these goals not just abstractly and 
combinatorially, but effectively, in a manner amenable to algorithmic computation. 
References such as |AA()S971 IA()92[ l(]H961 |Kap99t IMMP87I IMou85[ ESHi, which 
have their roots and applications in computational geometry, carry this out in the 
d = 2 case of boundaries of 3-polytopes (and for the first goal, on any polyhedral 
surface of dimension d = 2). Here, in arbitrary dimension d, our Theorem 15.21 savs 
precisely how past wavefront evolution determines the location in time and space of its 
next qualitative change. The combinatorial nature of Theorem 15 . 21 leads immediately 
to Algorithm 16.11 for effectively unfolding boundaries of polyhedra. 

The results and proofs in SectionsGHBlfor boundaries of convex polyhedra almost all 
hold verbatim in the more abstract setting of what we call (i-dimensional convex poly- 
hedral pseudomanifolds. The study of such spaces is suggested both by Stone's point 
of view in |Sto76j and by the more general methods in |BGP92j . Our Corollarv 17.121 
says that all convex polyhedral pseudomanifolds can be represented as quotients of 
Euclidean (usually nonconvex) polyhedral balls by identifying pairs of boundary com- 
ponents isometrically. The reader interested solely in this level of generality is urged 
to begin with Section [7| which gives a guide to Sections HHHl from that perspective, 
and provides the slight requisite modifications where necessary. Hence the reader can 
avoid checking the proofs in the earlier sections twice. 

The results in Section [7| on convex polyhedral pseudomanifolds are in many senses 
sharp, in that considering more general spaces would falsify certain conclusions. We 
substantiate this claim in Section |H1 where we also discuss extensions of our methods 
that are nonetheless possible. For example, we present an algorithm to construct 
geodesic Voronoi diagrams on boundaries of convex polyhedra in Section 18.91 

The methods of this paper suggest a number of fundamental open questions about 
the metric combinatorics of convex polyhedra in arbitrary dimension, and we present 
these in Section IHl Most of them concern the notion of vistal tree in Definition 19.11 
which encodes all of the combinatorial types of shortest paths (or equivalently, all 
bifurcations of the wavefront) emanating from a source point. The first two questions, 
Coniectures 19.21 and 19. 3t concern the complexity of our unfolding algorithm and the 
behavior of geodesies in boundaries of polyhedra. Along these lines, we remark also 
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on the complexity of nonconvex polyhedral manifolds, in Proposition 19.71 Our third 
question is about the canonical subdivision of the boundary of any convex polyhedron 
determined by the sets of source points having isomorphic vistal trees (Definition 1231 
and Conjecture 19 .5^ : it asks whether this vistal subdivision is polyhedral, and how 
many faces it has. Our final question asks how to realize unfoldings of polyhedral 
boundaries by embedded homotopies f Conjecture I9.11|) . 

Methods 

This section contains an extended overview of the paper, including background and 
somewhat informal descriptions of the geometric concepts involved. 

Unfolding polyhedra. While unfolding convex polytopes is easy |Ale48j , construct- 
ing a nonoverlapping unfolding is in fact a difficult task with a long history going back 
to Diirer in 1528 |0'Rnnj . When cuts are restricted to ridges (faces of dimension d — 1 
in a polyhedron of dimension d + 1), the existence of such unfoldings is open even for 
polytopes in p'ROnl [ZTe95] . It is known that nonconvex polyhedral surfaces need 
not admit such nonoverlapping unfoldings [BD EKMSJ. 

In this paper we consider unfoldings of a more general nature — when cuts are al- 
lowed to slice the interiors of facets. Such unfoldings are known, but only for three 
dimensional polytopes |SS86| IX092| ICH961 IAAOS97j . In fact, two different (although 
strongly related) unfoldings appear in these and other references in the literature: 
the Aleksandrov unfolding (also known as the star unfolding) |Ale481 [A.092j . and the 
source unfolding |VP71l l5S86j . Unfortunately, the construction of Aleksandrov un- 
foldings fails in principle in higher dimension f Section 18. 4p . As we mentioned earlier, 
we generalize the source unfolding construction to prove that the boundary S of any 
convex polyhedron of dimension d + 1, and more abstractly any convex polyhedral 
pseudomanifold S, has a nonoverlapping polyhedral unfolding U in W^. The second 
of the following two foldouts of the cube is a c? = 2 example of a source unfolding. For 




Figure 1. An edge-unfolding and a source unfolding of a cube into 

clarity, we present the discussion below in the context of boundaries of polyhedra. 

Cut loci. The idea of the source unfolding in arbitrary dimension d is unchanged 
from the case d = 2 of convex polyhedral surfaces. Pick a source point v interior to 
some facet (d-dimensional face) of S, so the tangent space is well-defined. Then, 
treating S like a Riemannian manifold, define the exponential map from Ty to S by 
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flowing along geodesies emanating from v. Our main unfolding result, Theorem 13.51 
says that exponentiation takes a certain open polyhedral ball C isometrically 
to a dense open subset of S consisting of points possessing a unique shortest path 
(length-minimizing geodesic) to v. The image of the closure of the open ball is 
all of S. The boundary \ maps onto the cut locus Ky, which by definition is the 
closure of the set of points in S with more than one shortest path to the source point v. 
These properties characterize Uy. 

In Riemannian geometry, when the manifold and the metric are both smooth, 
describing the cut locus for a source point is already an important and interesting 
problem (see |Kob89j for an excellent introduction and numerous references), although 
of course the exponential map can only be an isometry, even locally, if the metric is 
flat. Extending the notion of cut locus from Riemannian geometry to the polyhedral 
context is just as easy as extending it to arbitrary metric spaces. But showing that 
the open ball is a polyhedral foldout requires strong conditions on the complement 
of the cut locus, such as metric flatness and polyhedrality. We prove these results 
in Sections Q and |21 using methods based on the foundations of polyhedral geometry, 
and on Voronoi diagrams, culminating in Theorem 12.91 and Corollary 12.111 These 
conclusions depend crucially on convexity and do not hold in the nonconvex case. 

Geometry of wavefront expansion. Our existence proof for polyhedral nonover- 
lapping source foldouts, even given their Voronoi characterization in Theorem 12.91 
does not by itself provide a satisfactory combinatorial picture of the dynamics of 
wavefront expansion on polyhedra. For this, we must gain control over how the expo- 
nential map behaves as it interacts with warped points in S, namely those of infinite 
curvature, or equivalently points on faces of dimension d — 2 or less. 

Imagine the picture kinetically: the source point v emits a signal, whose wavefront 
proceeds as a (rf — l)-sphere of increasing radius — at least until the sphere hits the 
boundary of the facet containing v. At that stage, the wavefront folds over a ridge, 
or face of dimension d—1. Metrically, nothing has happened: points interior to 
ridges look to the wavefront just as flat as points interior to facets. But later, as the 
wavefront encounters faces of lower dimension, it is forced to bifurcate around warped 
points and interfere with itself, as signals emitted originally in different directions 
from V begin to curl around the infinite curvature and converge toward the cut locus. 

The question becomes: What discrete structure governs evolution of the wavefront 
on polyhedra? The most obvious first step is to define a finite collection of 'events', 
representing the points in time and space where the wavefront changes in some non- 
trivial way. If this is done properly, then it remains only to order the events according 
to the times at which they occur. However, in reality, the definition of an event is 
rather simple, while the geometry dictating time order of events is more complex. 

Starting from scratch, one might be tempted (and we were) to mark an event every 
time the wavefront encounters a new warped face. Indeed, this works in dimension 
d = 2 jlMEHIj: since the wavefront is a curve, its intersection with the set of edges 
is a finite set, and it is easy to detect when one of these points hits a vertex of S 
before another. But because the geometry is substantially more complicated in higher 
dimensions, in the end we found it more natural to say an event has occurred every 
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time the wavefront encounters a new facet through the relative interior of a ridge 
(see Definitions 12.31 and I4.14j) . This may seem counterintuitive, since the wavefront 
only interacts with and curls around faces of smaller dimension. However, wavefront 
collisions with warped points lead to intersections with ridge interiors infinitesimally 
afterward. In other words, the closest point {event point) on a facet to the source 
point V need not lie interior to a ridge, but can just as easily be warped. 

Again think kinetically: once the wavefront has hit a new face (of small dimension, 
say), it begins to creep up each of the ridges containing that face. Although in a 
macroscopic sense the wavefront hits all of these ridges simultaneously, it creeps up 
their interiors at varying rates. Therefore the wavefront hits some of these ridges 
before others in an infinitesimal sense. The moral is that if one wants to detect 
curling of the wavefront around warped faces, it is simpler to detect the wake of this 
interaction infinitesimally on the interiors of neighboring ridges. Sufficiently refined 
tangent data along ridges then discretizes the finite set of events, thereby producing 
the desired 'metric combinatorics' of wavefront expansion. 

Source poset. Making the above moral precise occupies Section |3J To single out a 
ridge whose interior is engulfed by the wavefront at a maximal rate (thereby making 
it closer to the source point) essentially is to find a ridge whose angle with the corre- 
sponding signal ray emitted from the source is minimal. When d = 2, this means that 
we do not simply observe two signals hitting vertices simultaneously, but we notice 
also the angles at which they hit the edges containing those vertices. The edge form- 
ing the smallest angle with its signal ray is the earlier event, infinitesimally beating 
out other potential events. (That each angle must be measured inside some ambient 
facet is just one of the subtleties that we gloss over for now.) 

To distinguish events in time macroscopically, only radii (distances from the source) 
are required. When d = 2, as we have just seen, a first derivative is enough to 
distinguish events infinitesimally. Generally, in dimension d > 2, one needs derivatives 
of order less than d, or more precisely, a directional derivative successively along each 
of d — 1 orthogonal directions inside a ridge. In Section ^ these derivatives are 
encoded not in single angles, but in angle sequences f Definition 14. 2p . which provide 
quantitative information about the goniometry of intersections between signal rays 
and the faces of varying dimension they encounter. More qualitative — and much more 
refined — data is carried by minimal jet frames fDefinition l4.1|) . which record not just 
the sizes of the angles, but their directions as well. 

The totality of the (finite amount of) radius and angle sequence data induces a 
partial order on events. The resulting source poset f Definition I4.14|) . which owes its 
existence to the finiteness result in Theorem 14.111 describes precisely which events 
occur before others — both macroscopically and infinitesimally. Since wavefront bi- 
furcation is a local phenomenon at an event point, incomparable events can occur 
simultaneously, or can be viewed as occurring in any desired order. Thus as time pro- 
gresses, wavefront expansion builds the source poset by adding one event at a time. 

The algorithm. It is one thing to order the set of events after having been given all 
of them; but it is quite another to predict the "next" event having been given only 
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past events. That the appropriate event to add can be detected locally, and without 
knowing future events, is the content of Theorem 15. 21 Its importance is augmented by 
it being the essential tool in making our algorithm for constructing the source poset, 
and hence also the source unfolding f Algorithm 16. Surprisingly, our geometric 
analysis of infinitesimal wavefront expansion in Sections |3] and El allows us to remove 
all calculus from Theorem 15.21 and hence Algorithm 16.11 detecting the next event 
requires only standard tools from linear algebra. 

As we mentioned earlier, our original motivation for this paper was its algorith- 
mic applications. Using the theoretical definitions and results in earlier sections, we 
present pseudocode (a semi-formal description) for our procedure constructing source 
unfoldings in Algorithm 16.11 That our algorithm provides an efficient method to 
compute source unfoldings is formalized in Theorem 16.51 

There are several reasons in favor of presenting pseudocode. First, it underscores 
the explicit effective nature of our combinatorial description of the source poset in 
Theorem 15.21 Second, it emphasizes the simplicity of the algorithm that results from 
the apparenly complicated analysis in Sections^^SJ in particular, the reader interested 
only in the computational aspects of this paper can start with Section Inland proceed 
backwards to read only those earlier parts of the paper addressed in the algorithm. 
Finally, the pseudocode makes Algorithm 16.11 amenable to actual implementation, 
which would be of interest but lies outside the scope of this work. 

A note on the exposition. Proofs of statements that may seem obvious based 
on intuition drawn from polyhedral surfaces, or even solids of dimension 3, demand 
surprising precision in the general case. Occasionally, the required adjustments in 
definitions and lemmas, and even in statements of theorems, were borne out only after 
considering configurations in dimension 5 or more. The definition of source image is 
a prime example, about which we remark in Section |H1 in the course of analyzing 
where various hypotheses (convexity, pseudomanifold, and so on) become essential. 
Fortunately, once the appropriate notions have been properly identified, the subtlety 
seemingly disappears: the definitions become transparent, and the proofs remain 
intuitive in low dimension. Whenever possible, we use figures to clarify the exposition. 

1. Geodesies in polyhedral boundaries 

In this paper, a convex polyhedron F of dimension c? is a finite intersection of closed 
half-spaces in some Euclidean space M'^, such that F that does not lie in a proper 
affine subspace of M*^. The polyhedron F need not be bounded, and comes with an 
induced Euclidean metric. Gluing a finite collection of convex polyhedra by given 
isometries on pairs of codimension 1 faces yields a (finite) polyhedral cell complex S. 
More precisely, S" is a regular cell complex endowed with a metric that is piecewise 
Euclidean, in which every face (closed cell) is isometric to a convex polyhedron. 

The case of primary interest is when the polyhedral cell complex S equals the 
boundary dP of a convex polyhedron P of dimension d + 1 in W^^^. 

Convention 1.1. We assume that S = dP is a polyhedral boundary in all theorems, 
proofs, and algorithms from here through Section |H1 
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We do not require P to be bounded, though the reader interested in polytopes will 
lose very little of the flavor by restricting to that case. Moreover, with the exception 
of Lemma 11.31 Proposition I2.1(J1 Corollary 12.111 and Theorem 13.51 the statements 
of all results from here through Section El are worded to hold verbatim for the more 
abstract class of convex polyhedral pseudomanifolds, as we shall see in Section [7| 

Denote by /i the metric on S, so /i(a, b) denotes the distance between points a,b & S. 
A path 7 C 5* with endpoints a and 6 is a shortest path if its length equals fi{a,b). 
Since we assume S has finitely many facets (maximal faces), such length-minimizing 
paths exist, and are piecewise linear. A path rj G S is a geodesic if i] is locally a 
shortest path; i.e. for every z & t] that is not an endpoint of rj, there exist points 
a, 6 G ?7 \ {z} such that 2; G 7 C for some shortest path 7 connecting a to b. 

Henceforth, as S has dimension d, a face of dimension d — 1 will be called a ridge. 
For convenience, we say that a point x is warped if x lies in the union 8^-2 of all faces 
in 5* of dimension at most d — 2, and call x flat otherwise. Every fiat point has a 
neighborhood isometric to an open subset of M.'^. 

Proposition 1.2. // a shortest path 'y in S connects two points not lying on a common 
facet, then 7 has no warped points in its relative interior. 

Proof. For any point w lying in the relative interior of 7, the intersection of 7 with 
some neighborhood of w consists of two line segments rj and rj' that are each straight 
with one endpoint at w, when viewed as paths in M.'^'^^. This is a consequence of local 
length-minimization and the fact that each facet of P is isometric to a polytope in M.'^. 
Moreover, if w happens to lie on a ridge while rj intersects the relative interior of some 
facet containing w, then local length-minimization implies that rj' is not contained in 
the facet containing 77. Lemma fl .31 shows that w does not lie in Sd-2, so the point w 
is not warped. □ 

Lemma 1.3. Let rj^rj' C S he two paths that (i) are straight in W^^^ , (ii) share a 
common warped endpoint w G Sd-2, o,nd (Hi) do not both lie in a single facet. There 
exists a neighborhood O of w in S such that for every a & t] (1 O and b E rj' r\0, the 
path rjab from a to w to b along rj and rj' is not a shortest path in S between a and b. 

Proof. Translate P so that w equals the origin G R'^'^^, and let Q be the unique 
minimal face of P that contains w. Since t] and 77' do not lie in a single facet, the 
2-plane E spanned by rj and rj' meets Q at exactly one point, namely 0. Since 
dim(Q) < d — 2, the span of Q and E has dimension at most d. Choose a line L 
whose direction is linearly independent from the span of Q and E. Then the 3-plane 
H = L + E intersects Q only at 0. Replacing P by P fl if, we can assume that 
dim(P) = 3, so that d = 2; note that is a vertex of if fl P by construction. 

Although the case d = 2 was proved in |Ale48l Theorem 4.3.5, p. 127] (see also 



jSS86l Lemma 4.1]), we provide a simple argument here, for completeness. Let O G S 
be the neighborhood of w consisting of all points at some fixed small distance from 
the vertex w. Then O can be laid flat on the plane by slicing along rj. One of the 
two points in this unfolding that glue to a G O connects by a straight segment in the 
unfolding to the unique point corresponding to b. This straight segment shortcuts rjab 
after gluing back to S. □ 
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Figure 2. Neighborhood of a vertex and its foldout after slicing along 
the segment rj. The points a and b are connected by a shortest path. 

Corollary 1.4. Let rj be a bounded geodesic in S starting at a point z not on any 
ridge. Then rj intersects each ridge in a discrete set, so rj traverses (in order) the 
interiors of a well-defined sequence Cri of facets (the facet sequence ofrj). 

Proof. Since rj is locally length minimizing, Proposition 11.21 implies that every inter- 
section of rj with a ridge takes place at a flat point. Such points have neighborhoods 
isometric to open subsets of M.'^, and these intersect rj in paths isometric to straight 
segments. It follows that rj intersects every ridge transversely. □ 

For each facet F of S = dP, let Tp be the affine span of F in M'^+^. 

Definition 1.5. Suppose two facets F and F' share a ridge R = F (1 F' . The folding 
map ^F,F' '■ Tp Tpi is the isometry that identifies the copy of R in Tp with the one 
in Tpi in such a way that the image of F does not intersect the interior of F' . 

In other words, the folding map $_f,f' is the rotation of Tp with [d— l)-dimensional 
axis R = F n F' so that F becomes coplanar with F' and lies on the other side of R 
from F'. It can be convenient to view as rotating all of R'^"'"^ instead of only 

rotating Tp onto Tp/. Informally, we say folds Tp along R to lie in the same 

affine hyperplane as F'. 

Definition 1.6. Given an ordered list C = [Fi, F2, . . . , Fi) of facets such that Fi 
shares a (unique) ridge with Fj+i whenever 1 < i < £, we write 

for the unfolding of Tp,, onto Tp-^, noting that indeed ^^^{TpJ = Tp^. Setting = 
(Fl, . . . , Fj), the sequential unfolding of a subset F C Fi U ■ ■ ■ U along C is the set 

(FnFi)U$^^^(FnF2)U---U<l>.i(FnF,) C Tp,, 
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By Corollary 11.41 we can sequentially unfold any geodesic. Next, we use this 
unfolding to show uniqueness of shortest paths traversing given facet sequences. 

Lemma 1.7. Let v and w be flat points in S. Given a sequence C of facets, there 
can be at most one shortest path 7 connecting v to w such that 7 traverses = C. 

Proof. Let 7 be a shortest path from f to w traversing C. Inside the union of facets 
appearing in C, the relative interior of 7 has a neighborhood isometric to an open 
subset of M'^ by Proposition 11.21 and the fact that the set of warped points is closed. 
Sequential unfolding of 7 into Tp for the first facet F in £ thus yields a straight 
segment in Tp. This identifies 7 uniquely as the path in S whose sequential folding 
along C is the straight segment in Tp connecting v to G Tp. □ 

In the proof of Lemma ITTTl we do not claim that the union of facets in the list C un- 
folds sequentially without overlapping, even though some shortest path 7 traverses C. 
However, some neighborhood of 7 in this union of facets unfolds without overlapping. 

Example 1.8. Consider the unfolding of a 1 x 1 x 3 rectangular box as in Figure El 
Denote by Fbot, -^top, -^front, -^back, -^left, -bright the bottom, top, front, back, left, and 
right faces, respectively. Denote by Ci the list of facets along which the points marked 
by i have been sequentially unfolded to create the foldout U C Tp^^^ in FigureEl Then: 

^1 = (-Pbot), ^2 = (-Fbot, -Fback), = (-Fbot , -Ffront ) , 

^4 = (-Pbot, -Pback, -Ptop), ^5 = (-Fbot , -Ffront , -Ptop) , ^6 = (-^bot , -Fback, -^left) , 

^7 = (-Fbot, -Ffront, -Pleft), = (-Fbot , -^back, -fright ) , ^9 = (-^bot , -Ffront , -fright ) • 




Figure 3. An unfolding of a 1 x 1 x 3 box. 



2. Cut loci 

Most of this paper concerns the set of shortest paths with one endpoint fixed. 

Definition 2.1. Fix a source point v ^ S lying interior to some facet. A point x & S 
is a cut point^ if x has more than one shortest path to v. Denote the set of cut points 
by K^, and call its closure the cut locus C S. 

^Ovlt: usage of the term 'cut locus' is standard in differential geometry, just as our usage of 'ridge' is 
standard in polyhedral geometry. However, these usages do not agree with terminology in computer 
science, such as in ^SS86. A092 : their 'ridge points' are what we call 'cut points'. Furthermore, 'cut 
points' in |A092| are what we would call 'points on shortest paths to warped points' (when d = 2). 
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Here is a consequence of Proposition 11.21 

Corollary 2.2. No shortest path in S to the source point v has a cut point in its 
relative interior. 

Proof. Suppose c is a cut point in the relative interior of a shortest path from v to w. 
Replacing the path from f to c with another shortest path from f to c yields a new 
shortest path from v to w. These two paths to w meet at the flat point c G S" by 
Proposition 11.21 The resulting Y-shaped intersection at c can be improved upon in a 
neighborhood of c isometric to an open set in (Fig. a contradiction. □ 




to V 



Figure 4. An intersection that is Y-shaped cannot locally minimize 
length in (segment 7 is a shortcut). 

Our study of polyhedrality of cut loci will use Voronoi diagrams applied to sets of 
points from the forthcoming definition, around which the rest of the paper revolves. 

Definition 2.3. Suppose that the source point v connects by a shortest path 7 to a 
point X that lies on a facet F or on one of its ridges R G F, but not on any face of 5* 
of dimension d — 2 or less. If the sequential unfolding of 7 into Tp is the segment [u, x], 
then u & Tp is called a source image for F. Let sicp be the set of source images for F. 

Lemma 2.4. The set sicp of source images for any facet F of S is finite. 

Proof. The shortest path in R'^^^ between any pair of distinct points x and y in a 
facet F is the straight segment [x, y] . Since this segment is actually contained in S, 
any shortest path 7 in S" must contain [x, y] whenever it contains both x and y. Taking 
X and y to be the first and last points of intersection between 7 and the facet F, we 
find that F can appear at most once in the facet sequence of a shortest path starting 
at the source point v. Hence there are only finitely many possible facet sequences of 
shortest paths in S. Now apply Lemma fl. 71 □ 

Example 2.5. Consider a unit cube with a source point in its bottom face, as in 
Fig El Then the top face has 12 source images, shown in FigEJ The four stars 'V 
are sequential unfoldings of the source point (along three ridges each) that are not 
source images: each point in the top face is closer to some source image than to any 
of these stars. 

The next result on the way to Theorem 12.91 generalizes |Mou85t Lemma 3.1] to 
arbitrary dimension. Its proof is complicated somewhat by the fact (overlooked in 
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Figure 5. Source point v on the 'bottom' face, 12 source images for 
the 'top' face of a cube and 4 'false' source images (view from the top). 



the proof of jMou85t Lemma 3.1]^) that straight segments can he inside the cut locus, 
and our lack of a priori knowledge that the cut locus is polyhedral. 

Proposition 2.6 (Generalized Mount's lemma). Let F be a facet of S, and suppose 
that u E STCp is a source image. If w E F, then the straight segment [i^, w] C Tp has 
length at least fi{v,w). 

Example 2.7. In Fig. IHl the left figure is a typical illustration of Proposition 12.61 in 





Figure 6. Generalized Mount's lemma (fails for the shaded region) 

dimension d = 2. In contrast, the right figure will never occur: any point w interior to 
the shaded region is closest to the source image u, but the straight segment connecting 
w to u has not been sequentially unfolded along the correct facet sequence. 

Proof. Since the two functions F —>■ W sending w to fi{v, w) and to the length of [z/, w] 
are continuous, we can restrict our attention to those points w lying in any dense 

^Much of |Mou85| . but not Lemma 3.1, was later incorporated and published in |MMP87] . 



METRIC COMBINATORICS OF CONVEX POLYHEDRA 



13 



subset of F. In particular, the cut locus has dense complement in F f Corollary 12. 2|) 
as does the boundary of F, so we assume throughout that w lies in neither the cut 
locus nor the boundary of F. 

Having fixed v G srcj?, choose a point x G F as in Definition 1231 The set srci?([a;, w]) 
of source images sequentially unfolded from shortest paths that end inside the segment 
[x, w] is finite by Lemma Hence we may furthermore assume that w does not lie on 
any hyperplane H that is equidistant from v and a source image v' G sicf{[x,w\). In 
other words, we assume w does not lie inside the hyperplane perpendicularly bisecting 
any segment [z/, z/'] . 

Claim 2.8. With these hypotheses, if v & sicp hut no shortest path unfolds sequen- 
tially to the segment [viw], then w is closer to some point v' G srcir([x,w]) than to v. 

Assuming this claim for the moment, we may replace v with v' and x with another 
point x' on [x, w\ . Repeating this process and again using that the set of source images 
sequentially unfolded from shortest paths ending in [x, w\ is finite, we eventually find 
that the unique source image uj G srcir([x, w]) closest to w is closer to w than v is. 
Since \jJ^w\ has length ij,{uj,w), it suffices to prove Claim 

Consider the straight segment [x, w], which is contained in F by convexity. Let Y 
be the set of points y G [x, w] having a shortest path •jy from v that sequentially 
unfolds to a segment in Tp with endpoint ly. Then Y is closed because any limit 
of shortest paths from v traversing a fixed facet sequence £ is a shortest path that 
sequentially unfolds along £ to a straight segment from the corresponding source 
image. Thus, going from x to w, there is a last point x' G Y. This point x' is by 
assumption not equal to w, so x' must be a cut point (possibly x = x'). 

There is a facet sequence C and a neighborhood O of x' in [x', w] such that every 
point in O connects to f by a shortest path traversing C, and such that unfolding 
the source along C yields a source image u' u in Tp. This point v' connects to x' 
by a segment of length /i(t>,x'), so the hyperplane H perpendicularly bisecting [i/, z/'] 
intersects [x, w] at x'. By hypothesis w ^ H, and it remains to show that w lies on 
the side of H closer to z/'. 

The shortest path from v to x' has a neighborhood in S disjoint from the set of 
warped points and hence isometric to an open subset of M"^ by Proposition ll.2t because 
x' is itself not a warped point (we assumed x lies interior to F or to a ridge R (Z F). 
After shrinking O if necessary, we can therefore ensure that each segment [u, y] for 
y G C is the sequential unfolding of a geodesic rjy in S. The geodesic rjy for y E 0\x' 
cannot be a shortest path by definition of x', so [z/, y] has length strictly greater than 
fi{v, y). We conclude that (9 \x', and hence also w, lies strictly closer to z/' than to u. 
This finishes the proof of Claim 12.81 and with it Proposition 12.61 □ 

Before stating the first main result of the paper, we recall the standard notion of 
Voronoi diagram V(T) for a closed discrete set T = {z/, z/', . . .} of points in M"^. This 
is the subdivision of R'' whose closed cells are the sets 

V{T, z/) = {C e M'^ I every point z/' G T satisfies |C - < |C - ^'\}- 

Thus ( lies in the interior of V{T, v) if is closer to v than to any other point in T. 
An example is given in Fig. [7| 
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Figure 7. An example of a Voronoi diagram. 

Theorem 2.9. Fix a facet F of S , and let Vd-i C Tp be the union of the closed 
cells of dimension d — 1 in the Voronoi diagram V(srci?) for the set of source images 
in Tp. If F° is the relative interior of F, then the set F° fl Ky of cut points in F° 
coincides with the intersection F° fl Vd-i- Moreover, if R° is the relative interior of 
a ridge R G F , then the set R° fl of cut points in R° coincides with R° fl Vd-i- 

Proof. Every shortest path from the source f to a point w in F° or R° unfolds to 
a straight segment in Tp of length ^{v,w) ending at a source image for F. Propo- 
sition 1221 therefore says that w lies in the Voronoi cell y(T,i/) if and only if the 
segment [z/, w] has length exactly fi{v,w). In particular, v has at least two shortest 
paths to w if and only if w lies in two such Voronoi cells — that is, w G V^-i. □ 

To illustrate Theorem 12.91 consider Example 12.51 The Voronoi diagram of source 
images gives the cut locus in the top face of the cube. 




Figure 8. Cut locus of the 'top' face of the cube. 

Theorem 12 . 91 characterizes the intersection of the cut locus with faces of dimension 
d or d — 1 in S. For faces of smaller dimension, we can make a blanket statement. 

Proposition 2.10. Every warped point lies in the cut locus K^; that is, Sd-2 ^ -^u- 
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Proof. It is enough to show that every point w in the relative interior of a warped 
face of dimension d — 2 is either a cut point or a hmit of cut points, because the cut 
locus is closed by definition. Let 7 be a shortest path from w to v. 

First assume that every neighborhood of w contains a point having no shortest 
path to V that is a deformation of 7. Suppose that (?/i)ieN is a sequence of such points 
approaching w, with shortest paths (7i)igN connecting the points i/i to v. Since there 
are only finitely many facets containing w and finitely many source images for each 
facet, we may assume (by choosing a subsequence if necessary) that for all i, the 
sequential unfolding of 7^ connects to the same source image for the same facet. The 
paths (7j)igN then converge to a shortest path 7' 7^ 7 to w from v, so w E K^. 

Now assume that every point in some neighborhood of w has a shortest path to v 
that is a deformation of 7. Every point on 7 other than w itself is fiat in S by 
Proposition E21 Therefore some neighborhood of 7 in 5 is isometric to an open subset 
of a product M''^^ x C, where C is a 2-dimensional surface that is fiat everywhere 
except at one point c G C (so C is the boundary of a right circular cone with apex c) . 
The set of points in M"'-^ X C having multiple geodesies to the image of f in M°' ^ x C 
is a relatively open half-space of dimension d — 1 whose boundary is M°'~^ x {c}. Some 
sequence in this open half-space converges to the image oi w. □ 

Theorem 12.91 and Proposition 12. lOl implv the following description of cut loci. 

Corollary 2.11. If v is a source point in S, then 

1. the cut locus Kf, is a polyhedral complex of dimension d — 1, and 

2. the cut locus is the union U Sci-2 of the cut points and warped points. 

Proof. Part |21 is a consequence of Theorem 12.91 and Proposition I2.1()| the latter tak- 
ing care of Sd-2, and the former showing that points in the cut locus but outside 
of Sd-2 are in fact cut points. Theorem 12.91 also implies that the intersection of the 
cut locus Ky with any closed facet F equals the polyhedral complex F fl Vd-i of di- 
mension d — 1. Using Proposition 12 . 1 01 in addition to the last sentence of Theorem 12. 91 
we conclude that these polyhedral complexes glue to form a polyhedral complex of 
dimension d — 1 that equals the cut locus K^, proving part[T] □ 



3. Polyhedral nonover lapping unfolding 

In this section we again abide by Convention 11.11 so S* is a the boundary of convex 
polyhedron P of dimension + 1 in W^^^. 

Definition 3.1. A polyhedral complex C S" of dimension — 1 is a cut set if K 

contains the union Sd-2 of all closed faces of dimension d — 2, and \ i^' is open and 
contractible. A polyhedral unfolding of 5* into M'^ is a choice of cut set K and a map 
S \ K ^M'^ that is an isometry locally on S' \ -ft". A nonoverlapping foldout of S is 
a surjective piecewise linear map ip : U S such that 

1. U is the closure of its interior U, which is an open topological ball in R'', and 

2. the restriction of (f to U is an isometry onto its image. 
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Note that K is not required to a polyhedral subcomplex of S, but only a subset 
that happens also to be a polyhedral complex; thus K can 'slice through interiors 
of facets'. The open ball U in item 1 of the definition is usually nonconvex. The 
polyhedron P is a polytope if and only if f/ is a closed ball — that is, bounded. 

When the domain ?7 of a nonoverlapping unfolding happens to be polyhedral, mean- 
ing that its boundary f/ \ f/ is a polyhedral complex, the image K = if{U \ U) is 
automatically a cut set in S. Indeed, piecewise linearity of (p implies that is a poly- 
hedral complex of dimension d — 1; while the isometry implies that K contains Sd-2, 
and that the open ball U = S \ K is contractible. Therefore: 

Lemma 3.2. // U is polyhedral, then a nonoverlapping foldout (f : U ^ S yields an 
ordinary polyhedral unfolding by taking the inverse of the restriction of ip to U . 

This renders unambiguous the term polyhedral nonoverlapping unfolding. 

The points in S outside of the [d — 2)-skeleton Sd-2 constitute a noncompact flat 
Riemannian manifold 5*°. When a point w lies relative interior to a facet F, the 
tangent space Tyj is identified with the tangent hyperplane Tp of F, but when w hes 
on a ridge, there is no canonical model for T^. 

Most tangent vectors Q ETyj can be exponentiated to get a point exp(C) G 5*° by the 
usual exponential map from the tangent space to the Riemannian manifold S° . 
(One can show that the set of tangent vectors that cannot be exponentiated has 
measure zero in T^; we shall not use this fact.) In the present case, we have a partial 
compactification S" of 5"°, which allows us to extend this exponential map slightly. 

Definition 3.3. Fix a point w & S° = S \ Sd-2- A tangent vector Q E Tyj can 
he exponentiated if the usual exponential of tC, exists in 5*° for all real numbers t 
satisfying < t < 1. In this case, set exp((^) = lim^^i exp(t^). 

The exponential map fc^ : t ^ exp(tC) takes the interval [0, 1] to a geodesic rj C S, 
and should be thought of as 'geodesic flow' away from w with tangent (. 

Henceforth fix a source point v E S not lying on any face of dimension less than d. 

Definition 3.4. The source interior consists of the tangent vectors ( eT^ at the 
source point v that can be exponentiated, and such that the exponentials exp{t() for 
< t < 1 do not lie in the cut locus K^. The closure of [7^, is the source foldout Uy. 

Our next main result justifies the terminology for and its closure Uy. 

Theorem 3.5. Fix a source point v in S . The exponential map exp : Uy ^ S 
from the source foldout to S is a polyhedral nonoverlapping foldout, and the boundary 
Uy \ Uy maps onto the cut locus Ky. Hence Ky is a cut set inducing a polyhedral 
nonoverlapping unfolding S \ Ky — > Uy to the source interior. 

Proof. It suffices to show the following, in view of parts ^ and |21 from Corollarv 12.111 

3. The metric space S \ Ky is homeomorphic to an open ball. 

4. The exponential map exp : Uy ^ S is piecewise linear and surjective. 

5. The exponential map exp : Uy ^ S \ Ky is an. isometry. 

Every shortest path is the exponential image of some ray in Uy by Proposition 11.21 
and the set of vectors ( E Uy mapping to S ^ Ky is star-shaped by part El along with 
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Proposition 11.21 and Corollary 12. 2[ This implies part El and surjectivity in part |3J 
The space 5*° = 5* \ Sd-2 is isometric to a fiat Riemannian manifold. Hence the 
exponential map is a local isometry on any open set of tangent vectors where it is 
defined. The definition of implies that exp is injective on the interior f4, so the 
surjectivity in part |3] shows that exp : f/^, — > 5 \ Ky is an isomorphism of Riemannian 
manifolds, proving part |S1 Every isometry between two open subsets of afiine spaces 
is linear, so the piecewise linearity in part 0] is a consequence of part El □ 

Example 3.6. Consider a cube P and a source point v located off-center on the 
bottom face of P, as in Example 12.51 and Figure |H1 The cut locus and the corre- 
sponding source foldout are shown in Fig. (HI See Fig. HI for the case when v is in 
the center of the bottom face. 



Remark 3.7. Surjectivity of the exponential map does not follow from 5"° being a 
Riemannian manifold: convexity plays a crucial role (see Fig.Elfor the case of a non- 



convex surface.) In fact, surjectivity of exp on a polyhedral manifold is equivalent — in 
any dimension — to the manifold having positive curvature |Sto76t Lemma 5.1]. The- 
orem 13.51 extends to the class of convex polyhedral pseudomanifolds, but not quite 
verbatim; see Theorem 17.111 for the few requisite modifications. 

4. The source poset 

In this section we define the source poset f Definition 14. 14|) . and in the next, we show 
how to build it step by step fTheorem l5.2|) . The reader should consider Definition 14. 141 




Figure 9. Cut locus K„ and source foldout U„ of the cube 




Figure 10. Shaded region lies outside of exp(T^). 
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as the main result in this section, ahhough it is the existence and finiteness properties 
for minimal jet frames in Theorem 14 . 1 1 1 1 hat endow the source poset with its power 
to make continuous wavefront expansion combinatorially tractable. 

Definition 4.1. Fix a polyhedron V in M°'. Given a list ( = {(i, . . . , (r) of mutually 
orthogonal unit vectors in Mf^, define for e G M the unit vector 

\i X & V and x + eJ(^{e) lies in V for all small e > 0, then the vector- valued function 
is a unit jet of order r at x in V , and C is a partial jet frame at x along V . If, in 
addition, x + eJ^{e) lies relative interior to V for all small e > 0, then is a, jet frame. 

The definition will be used later in the case where the convex polyhedron y is a 
closed Voronoi cell R fl VisiCF^uj) for some ridge i? of a facet F, and uj e src^ is 
a source image. Think of the point x & V a,s the closest point in V to uj. It will 
be important later (but for now may help in understanding the next definition) to 
note that the relative interior of a polyhedron V = Rr\ V{sicf, oj) is contained in the 
relative interior of the ridge R by Definition 12.31 and Theorem 12.91 

We do not assume the polyhedron V has dimension d. However, the order r of a 
unit jet in V , or equivalently the order of a jet frame along V , is bounded above by 
the dimension of V . In particular, we allow dim(\^) = 0, in which case the only jet 
frame is empty — that is, a list of length zero — and = 0. 

The lexicographic order on real vectors a and h of varying lengths is defined by 

(ai,...,a,.) < 

if the first nonzero coordinate of a — 6 is negative, where by convention we set = 
for i > r + 1 and hj = for j > s + 1. 

Definition 4.2. Fix a convex polyhedron V in M'^, a point x & V, and an outer 
support vector z/ G M"' for V at x, meaning that u ■ y < u ■ x for all points y & V. 
A jet frame ( aX x along V is minimal if the angle sequence — (z/ ■ Ci, . . . , z/ ■ (^r) is 
lexicographically smaller than — {u ■([,... ,1/ ■ C^,) for any jet frame (' at x along V. 

Again think of ^ = i? fl V{stcf, uj), with u = uj — x being the outer support vector. 

In general, that u is an outer support vector at x means equivalently that x is 
the closest point in ^ to a; + z/. Minimal jet frames ( can also be described more 
geometrically: the angle formed by u and (i must be as small as possible, and then 
the angle formed by z/ and ^2 niust be as small as possible given the angle formed by z/ 
and (i, and so on. It is worth bearing in mind that because z/ is an outer support 
vector, the angle formed by z/ and (i is at least tt/2 (that is, obtuse or right). 

Lemma 4.3. // ( and (' are vectors of equal length in M'^, and z/ G M'^ a vector 
satisfying i' ■ ( < and ly ■ (' < 0, then — (\ < — ('\ if and only if u ■ ( > u ■ (' . 

Proof. Draw z/ pointing away from the center of the circle containing ( and with 
these vectors on the other side of the diameter perpendicular to u. Then use the 
law of cosines: the radii ( and (' have equal length, and u has fixed length; only the 
distances from ly to ( and (' change with the angles of ( and (' with z/ (see Fig.lll|). □ 
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Figure 11. Dot product vs. length 

Minimal jet frames admit a useful metric characterization as follows. 

Proposition 4.4. Fix two polyhedra V and V with outer support vectors u and v' , 
of equal length, at points x E V and x' G V , respectively. Let ( and (' be partial 
jet frames at x along V and x' along V , respectively. The angle sequence —v ■ ( is 
smaller than —v' ■ (' in lexicographic order if and only if there exists Eq > such that 
X + V is closer to x + eJ^{e) than x + u' is to x' + eJ^i{e) for all positive e < sq. 

Proof. Since the dot product of z/ with each vector J^{e) or J^/(e) is negative, and these 
are unit vectors, it is enough by Lemma f4. 31 to show that minimality is equivalent to 

^ ■ J({^) > ^ ■ <^c'(^) for all nonnegative values of e < Eq. 

If the first nonzero entry of u ■ ( — u ■ (' is c = i^-id — Q), then for nonnegative values 
of e approaching 0, the difference u ■ J^{b) — v ■ Jqi{s) equals ce*~^ times a positive 
function approaching 1. The desired result follows easily. □ 

Corollary 4.5. Fix an outer support vector v at a point x in a polyhedron V . A jet 
frame Q at x along V is minimal if and only if, for every jet frame (' at x along V, 
X + u is weakly closer to x + eJ^{e) than to x + eJ^,{e) for all small nonnegative e. 

It is not immediately clear from the definition that minimal jet frames always 
exist: a priori there could be a continuum of choices for Ci, and then a continuum of 
choices for ^2 in such a way that no minimum is attained. Although such continua of 
choices can indeed occur, we shall see by constructing minimal jet frames explicitly 
in Theorem 14.111 that a minimum is always attained. 

First we need to know more about how (partial) jet frames at x reflect the local 
geometry of V near x. The tangent cone to a polyhedron V C M'^ at x G is the cone 

%V = RyoiCeW'lx + CeV} 

generated by vectors that land inside V when added to x. 

Definition 4.6. Fix a partial jet frame ^ at x along a polyhedron V in W^. Let (-^ be 
the linear subspace of M"^ orthogonal to the vectors in (, and fix a sufficiently small 
positive real number e. Then define the iterated tangent cone 

T^V = T^{{^ + C^)nT,V) 
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as the tangent cone at ,^ = J(^{e) to the intersection of T^V with the affine space ^ + C"*"- 

Just as the partial jet frames of order 1 generate the tangent cone T^V, we have 
the following characterization of iterated tangent cones. We omit the easy proof. 

Lemma 4.7. The iterated tangent cone T^V is generated by all unit vectors (r+i in 
extending the partial jet frame ( = (Ci, . . . , Cr) to a partial jet frame (Ci, ■ ■ ■ ,Cr, C+i) 
of order r + 1. In particular, iterated tangent cones do not depend on the small e > 0. 

Now we set out to construct minimal jet frames inductively. 

Lemma 4.8. Fix a polyhedron V and an outer support vector u at x & V . If C ^ T^V 

is a unit vector with v ■ ( maximal, then v is an outer support vector at & T^V. 

Proof. If C £ T^'^ is a unit vector satisfying u ■ > 0, then ^ = (C + ^C) I VT+T^ for 
small £ > is a unit vector in satisfying v-^ > u-^, contradicting maximality. □ 

In 'generic' cases, the functional C ' C fo^^ outer support vector z/ on a cone 
will take on the maximum value zero uniquely at the origin. In this case, as we now 
show, there can be only finitely many unit vectors ( in the cone having u ■ ( maximal, 
and these lie along the rays, meaning one- dimensional faces of the cone. Note that 
genericity forces the cone to be sharp, meaning that it contains no linear subspaces. 

Proposition 4.9. Let v he an outer support vector for a sharp polyhedral cone C , 
and assume v is maximized uniquely at the origin 0. The minimum angle between v 
and a unit vector occurs when ( lies on a ray of C. 

Proof. Let Z be the set of unit vectors in C. Suppose that L is a 2-dimensional 
subspace inside the span of C, and let z/ be the orthogonal projection of z/ onto L. 
View 1/ and z/ as functionals on L via C ' ^ind observe that u ■ ( = u ■ ( for all 
^ e L. The circular arc Z (1 L lies inside the unit circle in L, and u takes nonpositive 
values on Z n L because u is an outer support vector. Elementary geometry shows 
that z/ is therefore maximized on Z (1 L only at one or both of the endpoints of 
the arc Z n L. This argument proves that u cannot be maximized on Z at a point 
C ^ Z unless ( lies in the boundary of Z. The result now follows by induction on the 
dimension of the cone C. □ 

In 'nongeneric' cases, including when the polyhedral cone C has nonzero lineality, 
which is by definition the largest vector space contained in C, the functional u is 
maximized along a face of positive dimension. In this case, there is always a continuum 
of choices for unit vectors C & C having u ■ ( = 0. However, the sequences of iterated 
tangent cones to appear in Theorem 14. 1 II will not in any noticeable way depend on 
the continuum of choices, because of the next result. 

Lemma 4.10. Fix a polyhedron V, a point x E V, and a face F of V containing x. 
The iterated tangent cone T^V is independent of the jet frame ( for F at x. 

Proof. Translate V so x + 6J^{e) lies at the origin e R'^. Then F spans a dimension 
dim(F) linear subspace (F) C R'', and the iterated tangent cone is T^V = {F)-^nTQV. 
Now use the fact that TqV = T^V for all vectors ^ relative interior to F. □ 
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The main theorem in this section says that given an outer support vector u, there 
is a finite procedure using elementary hnear algebra for producing a single jet frame 
that is, in a precise sense, tilted as much toward v as possible. 

Theorem 4.11. Fix a polyhedron V and an outer support vector v at x (z V . In- 
ductively construct a finite set of jet frames for V at x by iterating the following 
procedure. For each of the finitely many partial jet frames ( already constructed: 

• If u is orthogonal to a nonzero vector in T^V , then add any such vector to C,. 

• If u ■ ( < for all nonzero vectors ( in T^V , then create one new partial jet 
frame for each of the (finitely many) rays ofT^V minimizing the angle with v, 
by appending to ( the unit vector along that ray. 

At least one of the finitely many jet frames constructed in this way is minimal. 

Proof. The sequences of vectors constructed by the iterated procedure are jet frames 
by Lemma f4. 81 Given an arbitrary jet frame ^ for ^ at x, it is enough to show that 
the angle sequence of ^ satisfies i/-^ > z/-^ in lexicographic order for some constructed 
jet frame (. Indeed, then a jet frame whose angle sequence is lexicographically min- 
imal among the constructed ones is minimal. Suppose that the first i — 1 entries 
(^1, . . . , agree with a constructed jet frame, but that (^i, . . . do not. 

If z/ ■ < then z/ ■ is less than z/ ■ Q for some constructed jet frame ( agreeing 
with ^ through the {i — entry, by Proposition 14.91 

If, on the other hand, z/ ■ = 0, then pick the index j maximal among those 
satisfying ^ and also z/ ■ = ■ ■ ■ = z/ ■ = 0. If there is a constructed jet frame C 
that agrees with ^ through the j**^ entry, but has v ■ ^j+i < v ■ Cj+i = 0, then already 
we are done. Therefore we can assume that the constructed jet frame C, agrees with ^ 
through index (z — 1), that has z/ ■ (^j = ■ ■ ■ = z/ ■ = 0, and that either v ■ Cj+i < 
or else C has order j. Replacing the vectors ^i, ■ ■ ■ in ^ with Q, . . . Xj yields a new 
jet frame by Lemma l4 . 1 01 applied to the face F of the iterated tangent cone T£V 
orthogonal to u, where (' = {(i, . . . ,Ci-i)- Downward induction on the number of 
entries of ^' shared with a constructed jet frame completes the proof. □ 

Our goal is to apply jets to define a poset structure on the set of source images. 
First, we need some terminology and preliminary concepts. The next definition is 
made in slightly more generality than required for dealing only with complete sets of 
source images because we shall need it for Theorem 15.21 

Resume the notation from previous sections regarding the polyhedral complex S. 
Recall that Tp = M'^ is the tangent hyperplane to the facet F. Removing from Tp 
the affine span Tr of any ridge R <Z F leaves two connected components (open half- 
spaces). Thus it makes sense to say that a point v (iTp lies either on the same 
side or on the opposite side of R as does F. 

Definition 4.12. Fix a facet F, a ridge R G F, and a finite set T C Tp. 

1. A point u; G T can see F through R in V(T) if uj lies on the opposite side of R 
as F does, and the closed Voronoi cell V{T,uj) contains a point interior to R. 

2. A point G T can see R through F in V{T) ii u lies on the same side of R as 
F does, and the closed Voronoi cell V{T,u) contains a point interior to R. 
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3. In either of the above two cases, the ridge R lies at radius r = r{R,uj) from uj 
if r equals the smallest distance in Tp from to a point of Rn V{T, u). 

4. The unique closest point p{R,uj) to uj in RCi V{T,uj) has distance r from u. 

5. The outer support vector of the pair {R,uj) is u — p{R,uj). 

6. The angle sequence Z{R,u) is the angle sequence —{u — p{R,uj)) ■ C for any 
minimal jet frame ( at p{R, uj) along R fl V{T, uj). 

Example 4.13. Fig. depicts examples of the notions from Definition 14.121 The 



solid pentagon is the face F, while the set T contains four points. The point on can see 
F through the ridge R, and can see the ridges R' as well as R" through F. The three 
closest points for these are indicated, as is the outer support vector for {R,uj). The 
point u' can see the ridge R' through F, but lj' can not see R" through F, because lu 
is closer to every point of R". 

In our applications, the finite set T will always be a subset of source images in srcp, 
often a proper subset. Now we are ready for the main definition of this section. It may 
help to recall that each source image u G srci;' can see F through a unique ridge R 
by Theorem 12.91 when M.'^ = Tp and the finite set T in Definition 14. 121 equals srci;'. 

Definition 4.14. Fix a source point v in S. An event is a pair (z/, F) with v G sicp 
a source image for the facet F . The event (z/, F) has 

1. radius r(i/, F) equal to the radius r(i?, v) from v to the ridge R through which 
V can see F in the Voronoi subdivision V(srci?) of Tp\ 

2. event point p{y, F) equal to the closest point z/) in Rf\V {sicp, z/) to z/; and 

3. angle sequence ^{y^F) equal to the angle sequence Z(i?, z/). 

(The trivial event {v, facet (f)) has radius 0, event point and empty angle sequence.) 
The source poset src(f , S) is the set of events, partially ordered with (z/, F) -< (z/', F') if 

• r{v,F) < r(z/',F'), or if 

• r{v,F) = r{v',F') and Z(z/, F) is lexicographically smaller than Z(z/',F'). 




Figure 12. Illustrations for Definition 14.121 
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Remark 4.15. Corollary 14.51 savs that breaking ties by lexicographically comparing 
angle sequences at event points is the same as breaking ties by comparing distances 
from each source image to a minimal jet at its event point. This is the precise sense 
in which the source poset orders events by comparing infinitesimal expansion of the 
wavefront along the interiors of ridges containing event points. 

5. Constructing source images 

Aside from its abstract dynamical interpretation, the importance of the source poset 
here stems from its ability to be computed algorithmically, as we shall see here and in 
Sectional Source images are built one by one, using only previously built source im- 
ages as stepping stones. These stepping stones form an order ideal in src(f , S), mean- 
ing a subset X C src(f , S) closed under going down: E ^ I and E' -< E ^ E' E I. 

To make a precise statement in the main result, Theorem l5.2l we need one more dose 
of terminology, describing constructions in S determined by a choice of order ideal. 

Definition 5.1. Fix an order ideal X in the source poset src(v, S). For each facet F, 
let Tp C Tp be the set of source images uj G sicp with {uj,F) e X. The set Sj of 
potential events consists of triples {u, F, B!) such that 

• uj can see the ridge R' through F in the Voronoi diagram V(T^), but 

• a second facet F' contains R' , and the unfolding u' = ^p^p/{uj) of u onto the 
tangent space Tp/ results in a pair (cu', F') that does not lie in X. 

If {u', F') is an event in src(f , S) \X, then we say it is obtained by processing {uo, F, R'). 
A potential event E E Sj is minimal if it has minimal radius r among potential events, 
and lexicographically minimal angle sequence among potential events with radius r. 

Tracing back through notation, if E = {u, F, R') is a minimal potential event, then 
the minimal radius is r = r(i?', z/), and the minimal angle sequence is Z{R', v). 

Theorem 5.2. Given a nonempty order ideal I in the source poset sic{v,S), pick a 
minimal potential event {v,F^R') in Ex- If i^' = ^p^p>{i') is the unfolding of u to the 
other facet F' containing R' , then X' = X U {{i^', F')} is an order ideal in src(f , S). 

The statement has two parts, really: first, u' G Tpi is indeed a source image; and 
second, X' is an order ideal in the poset src(f , S). To prove the theorem we need a 
number of preliminaries. We state results requiring an order ideal inside the source 
poset src(f , S) using language that assumes an order ideal X has been fixed. 

Recall from Section^the notion of facet sequence for a shortest path 7. If, on the 
way to a facet F', a shortest path 7 from the source point v traverses a facet F, then 
the corresponding source images in F and F' have a special relationship. Precisely: 

Definition 5.3. Let (z/, F) -< (i/', F') be events in the source poset. Suppose some 
shortest path 7 has facet sequence = (Fi, . . . , F^/) with a consecutive subsequence 

C = {Fe,..., Fe) in which F = F^ and F' = F^. 

If u' = $£(z/) = $£^(f) is the sequential unfolding of the source along 7, and also the 
sequential unfolding oi u E Tp into Tpi, then (i^, F) geodesically precedes (z/',F'). We 
also say that the shortest path 7 described above is geodesically preceded by {i^,F). 
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Since the Voronoi cells in Theorem 12.91 come up so often, it will be convenient to 
have easy terminology and notation for them. 

Definition 5.4. Given a source image uj G srci?, the cut cell of is = V{sicfiOj). 

Roughly speaking, our next result says that angle sequences increase at successive 
events along shortest paths, when the event point is pinned at a fixed point x. 

Proposition 5.5. If{v^F) geodesically precedes (z/',F') then (z^, -F) -< {v',F'). 

Proof. Because of the way partial order on src(f , S) is defined, we may as well assume 
that F and F' share a ridge i?', and that v' = $i?^i7"(z/) is obtained by folding along this 
ridge. In addition, we may as well assume that both event points p(z/, F) and p(z/', F') 
equal the same point x G S*, since otherwise r{i>,F) < r{i>',F'). Translate to assume 
this point x equals the origin 0, to simplify notation. Let R be the ridge through 
which z/ can see F, and set V = RdVu and V = R' H Vyi] these are the cut cells 
through which the source images v and v' see their corresponding facets. 

The angle geometry of u' relative to V in Tpi is exactly the same as the geometry 
of V relative to V in Tp, because v' is obtained by rotation around an axis in W^^^ 
containing V . In other words, v — v' is orthogonal to V . Therefore we need only 
compare the angles with u of jets along V and V . All jet frames will be at x. 

Suppose the finite sequence (^i, ^2, . . .) is a jet frame along V . Noting that V and 
V have disjoint interiors, choose the index r so that ^ = (^1, . . . ,^r-i) is a partial 
jet frame along V, but ^' = is not. It is enough to demonstrate that 

some partial jet frame (^1, . . . ,^r-i,Cr) along V has lexicographically smaller angle 
sequence than S,' ■ Equivalently, it is enough to produce a unit vector Cr in the iterated 
tangent cone T^V satisfying v ■ Cr > ■ ^r- 

Since R' CiVu' = R' (1 Vy by Theorem 12. 91 every line segment from z/ to a point in V' 
passes through V . Therefore, since we have translated to make a; = 0, every segment 
connecting v to T.J.V' passes through T^^V . This observation will become crucial below; 
for now, note the resulting inequality dim(l^) > dim(\^'), which implies that the 
iterated tangent cone T^V contains nonzero vectors. All such vectors by definition lie 
in the subspace orthogonal to the space (^) with basis ^1, . . . , ^r-i- The same holds 
for so we may replace v with a vector G by adding a vector in (^), since then 

■ C = ^ ' C for all vectors C, G S,'^. 

Fix a small positive real number e. The line segment [z/, J^i{sy\ intersects T^V at 
a point near J^{s). The image segment in by orthogonal projection modulo (^) is 
[cj, A^r], for A = e''/ ^/e"^ + ■ ■ ■ + e"^^. This image segment passes through the cone T^V 
at some point C on its way from uj to A^^- Elementary geometry of the triangle with 
vertices 0, u, and A^^ (see Fig. IT^ shows that the angle between u and ( is smaller 
than the angle between uj and X^r- Taking (r = C/ICI completes the proof. □ 

After choosing a minimal potential event E, we must make sure that when all is 
said and done, none of the other potential events end up below E in the source poset. 

Lemma 5.6. Suppose {uj, F, R') & Sj is a potential event with angle sequence Z and 
radius r. Let F' he the other facet containing R' and uj' = ^p^pi{uj) the unfolding of 
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Figure 13. Geodesic precedence implies smaller angle sequence 

uj ETp onto Tpi. If {uj', F') is an actual event, then it either has radius strictly bigger 
than r, or else its angle sequence Z{uj',F') is lexicographically larger than Z. 

Proof. Assume that {u', F') is an event. Quite simply, the result is a consequence of 
the fact that the cut cell -R' fl K;/ = i?' fl must be contained inside R! fl V{Tf, cu), 
which follows because Tp C sicp- □ 

In comparing a newly processed event (in source poset order) to other as yet un- 
processed events, we will need to know approximately how those other events will 
eventually arise. This requires the forthcoming lemma, in which a flat triangle in- 
side 5* is any subset of 5* isometric to a triangle in the Euclidean plane M^. 

Lemma 5.7. Fix a point x & S. There is an open neighborhood Ox of x in S such 
that, given y G Ox and a shortest path 7 from the source point v to y, some shortest 
path 7' from v to x has the following property: the loop formed by traversing 7' and 
then the segment [x,y] and finally the reverse of bounds aflat triangle in S. 

Proof. Choose Ox so small that the only closed faces of the cut locus intersect- 
ing Ox are those containing x. Every cut cell containing y & Ox also contains x by 
construction. Convexity of cut cells (Theorem 12. 9|) implies that the segment [x, y] 
lies inside every cut cell containing y (there may be more than one if y is itself a cut 
point). The source image obtained by sequentially unfolding 7 therefore connects to 
every point of [x, y] by a straight segment that sequentially folds to a shortest path. 
The union of these shortest paths is the flat triangle in question. □ 

Conveniently, all of the shortest paths to x already yield events in X: 

Lemma 5.8. Suppose some minimal potential event E E £j has closest point x. Let 
G be the last facet whose interior is traversed by a shortest path 7 from the source point 
to X. If UJ E stcg is the source image sequentially unfolded along 7, then {uj,G) G X. 

Proof. As 7 enters G, it crosses the relative interior of some ridge of G at a point w. 
The event point p{ij, G) can be no farther than w from u. On the other hand, 
fj,{v,w) < fj,{v,x), because 7 traverses the interior of G. Therefore {uj,G) has radius 
less than r{E) = fi{v,x). □ 

Proof of Theorem 15. 2L Suppose the minimal potential event {u, F, R') has closest 
point X = p{R', z/) to R' n V{Tf, z^), of radius r, and a minimal jet frame ( a.t x with 
angle sequence Z. 
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Let 7 be a shortest path from the source that ends at a point in the neighborhood 
from Lemma l^n By that lemma and Lemma 7 unfolds to produce a source image 
whose event either lies in X, or is obtained by processing a potential event in Sj, or 
is geodesically preceded by such a processed event. Applying Lemma 1^751 and then 
Proposition I5.5| we find that all events with event point x that are not in X have 
angle sequences lexicographically larger than Z. 

For positive e, set y{e) = x + eJ^{e). When e is small enough, y{e) lies interior 
to R', and close to x, in the neighborhood from Lemma 15 .71 By the previous 
paragraph, every source image containing y in its cut cell is either in X, or has angle 
sequence lexicographically larger than Z. 

Let us now compare, for all small positive e, the distance to y{e) from u with the 
distance to y{e) from any source image in srcp or srcp'. Clearly the distance from 
a source image u is minimized when y{e) lies in the cut cell V[j. Moreover, we may 
restrict our attention to those source images u whose cut cells Kj contain y{e) for all 
sufficiently small positive e. Definition 14.11 savs that C is a jet frame at x along V^j. 
Therefore, by Proposition 14.41 we conclude using the last sentence of the previous 
paragraph that y{6) is weakly closer to u than to u for all small positive e. This 
argument shows that z/' is a source image, so {u', F') is an event. Moreover, it shows: 

Claim 5.9. Any minimal jet frame ( at the event point x = p{R', v) along the poly- 
hedron R! n V{Tf, v) is a minimal jet frame at x along R' fl Vyi . 

Every event in src(f , 5*) \X is either obtained by processing a potential event in 
or is geodesically preceded by such a processed potential event. Using Claim we 
conclude by Lemma 15.61 and Proposition 15.51 that X' is an order ideal. □ 

6. Algorithm for source unfolding 

The primary application of the analysis up to this point is an algorithmic construction 
of nonoverlapping unfoldings of convex polyhedra, which we present in pseudocode 
followed by bounds on its running time. In particular, we show that the algorithm 
is polynomial in the number of source images, when the dimension d is fixed. (Later 
we shall state Conjecture 19.21 which posits that the number of source images is poly- 
nomial in the number of facets.) Other applications, some of which are further dis- 
cussed in SectionlSl include the discrete geodesic problem fCorollarv l6.6p and geodesic 
Voronoi diagrams f Algorithm 18. II in Section 

Roughly, Algorithm inn] consists of a single loop that with every iteration constructs 
one new event. Each event is a pair consisting of a facet and a point that we have called 
a source image in the affine span of that facet. The loop is repeated exhaustively until 
all of the events are computed, so the affine span of every facet has its full complement 
of source images. The Voronoi diagram for the set of source images in each affine 
span induces a subdivision of the corresponding facet. For each maximal cell in this 
subdivision, the algorithm computes a Euclidean motion (composition of rotation 
and parallel translation) that moves it into the affine span of the facet containing 
the source point. The union of these moved images of Voronoi cells is the output 
foldout in the tangent space T„ to the source point v. 
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At each iteration of the loop, the algorithm must choose from a number of potential 
events that it could process into an actual event. Each potential event E consists 
of an already-computed event {u, F) plus a ridge R in the facet F. Processing the 
event E applies a rotation to move the source image u into the affine span of the 
other facet containing R. The potential event that gets chosen must lie as close to 
the source point as possible; this distance is the radius r = r{R, u) at the beginning of 
the loop. The loop then calls Routine 16.21 to choose which event to process; although 
this routine is quite simple in structure, it is the part of the algorithm that most 
directly encounters the subtlety of working in higher dimensions. The end of the loop 
consists of updating the sets of source images and potential events; the latter requires 
Routine 16.31 which we have isolated because it is the only time-consuming part of the 
algorithm, due to its Voronoi computation. 

Let us emphasize that once a source point v is computed, it is never removed. This 
claim is part of Theorem 15.21 in which the correctness of Algorithm 16.11 — and indeed 
the procedure of the algorithm itself — is more or less already implicit, as we shall see 
in the proof of Theorem 16.41 

We assume that the convex polyhedron P is presented in the input of the algo- 
rithm as an intersection of closed half-spaces. Within the algorithm, we omit the 
descriptions of standard geometric and linear algebraic operations, for which we refer 
to |GU971 rPSSSj . These operations include the determination of lower-dimensional 
faces (such as ridges) given the facets of P, and the computation of Voronoi diagrams. 

Some additional notation will simplify our presentation of the algorithm. Denote 
by and IZ the sets of facets and ridges of P, respectively. If a ridge R E 71 lies in 
a facet F E J^, denote by 4>{F, R) the other facet containing P, so P fl 0(P, P) = P. 
Finally, for each facet F E J-', denote by Sp the set of all triples {u, P, P) such that u 
lies in the affine span Tp of P, and P G 7^ is a ridge contained in P. 

Algorithm 6.1. (Computing source unfolding) 

INPUT convex polyhedron P C M''"'"^ of dimension d + 1 

point V lying in the relative interior of a facet P of P 
OUTPUT source foldout of the boundary S = dP into = (see Sec. Ej) 
DEFINE for each P G JF: a finite set Tp C Tp of points 

for each pair (i^, P) satisfying u G Tp: an ordered list Ci^^p of facets 
for each P G JF: a set Sp C Sp of potential events 
£ = [Jp^-fSp, the set of all potential events 
INITIALIZE for P G JF: if w G P then Tp := and Sp = 0; 

otherwise := {v}, C^^p := (P), £p := {iv,F,R) | P G 7^ and P C P} 
COMPUTE ^p^pf for all F,F' e such that P n P' G 7?. is a ridge (see Def. [T3|l 

WHILE £ 

DO r := min{r(P, u) \ {u, P, P) G 8} (see Def. 

CHOOSE A POTENTIAL EVENT E = (z/, F, R) E £ TO PROCESS 
set P':=0(P,P), z/' := $^,^,(z/), C,,,p> := iC,,p F') 
update Tp/ ^ Tp/ U {u'} 

£p/ ^ {{uj,F',R') G £p/ such that uj G Tp/, and 
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POINT UJ e Tp/ CAN SEE R' THROUGH F' , and 

Lj' i Tg, where G = <i){F\R'), J = ^f',gH} 
£f *— £f ^ {E}, £ ^ {JaeT^G 

END WHILE-DO 

COMPUTE for all facets F e and points u G Tp: 

$£ for C = Cy^F (see Def. [LSI), then 
U^{v, F) := [F n ViTF, v)) C T„ (see Thm. HHJ 
RETURN the foldout t4 = [J Uv{i^-, F), the union being over all F G and u E Tf 

Routine 6.2. (choose a potential event to process) 

INPUT the set £ = Ufc:^^^ °f potential events, and the radius r > 
OUTPUT an event E e £ (see Def. 15. ip 
COMPUTE the closest potential events £o := {{uj, F,R)e£\ r{uj, F) = r} 
angle sequence Z{R,uj) for all {uj,F,R) G £o (see Def. 14.12(1 
FIND a potential event E = (uj,F,R) G £o with lexicographically 
minimal angle sequence Z{R,uj) (see Sec. Ej) 
RETURN the event E = (cu, F, R) 

Routine 6.3. (point lu e T can see R through F) 

INPUT facet F G JF, ridge R E TZ, finite set of points T C Tp, and a; G T 
OUTPUT boolean variable 3 G {True, False} (see Def. I4.12p 
COMPUTE Voronoi diagram V(T) (see Sec. \^ 

IF Voronoi cell V{T,to) C V(T) contains a point interior to R 
and CO lies on the same side of i? as F does in Tp 
then D := True; 
otherwise H := False 
RETURN the variable 3 

In the pseudocode, we have used the two different symbols and ':=' to dis- 
tinguish between those variables that are being updated and those that are being 
completely redefined at each iteration of the while-do loop. We hope this clarifies 
the structure of Algorithm 16.11 

Theorem 6.4. For every convex polyhedron P C R"'^-'^ with boundary S = dP , and 
any source point v in a facet of S, Alaorithm \6. 1\ computes the source foldout '^Ty. 

Proof. First, we claim by induction that after each iteration of the while-do loop, 
the set {(z/, F) I F G and z/ G Tf} is an order ideal in the source poset src(f , S) from 
Definition l4.14l The claim is clear at the beginning of the algorithm. By construction. 
Routine 16.21 picks a minimal potential event E to process. The loop then adds an 
event by processing F, with the aid of Routine 16.31 Theorem 15.21 implies that what 
results after processing F is still an order ideal of events, proving our claim. Since 
the poset src(f , S) is finite by Lemma the algorithm halts after a finite number of 
loop iterations. Finally, by Theorem 12.91 the Voronoi cells in each facet coincide with 
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the polyhedral subdivision of each facet by the cut locus K^^ so Theorem 13.51 shows 
that the foldout in the output is the desired (nonoverlapping) source foldout U^. □ 

For purposes of complexity, we assume throughout this paper that the dimension d 
is fixed. Thus, if the convex polyhedron P C M'^+^ of dimension d has n facets, so 
P is presented as an intersection of n closed half-spaces, we can compute all of the 
vertices and ridges of P in polynomial time |GU97t IZie95j . For simplicity, we assume 
these are precomputed and appended to the input. 

The timing of Algorithm 16 . II cruciallv depends on the number of source images. Let 

IfCy := max Isrci;' 

be the largest number of source images in a tangent plane Tp for a facet F. (This 
number can change if the source point v is moved. For example, src^ = 4 if f is in the 
center of a face, while srcy = 12 if t> is off-center as in Figures |H1 and El) Note that com- 
puting Voronoi diagrams for points in M'^ can be done in iV°('^) time |G097l p.381]. 
See |Aur9H ICha9H IFor95j for details and further references on Voronoi diagrams, 
and |(7()97[ IPS85j for other geometric and linear algebraic computations we use. 

Theorem 6.5. When the dimension d is fixed, the cost of Alaorithm \6. 1\ is polynomial 
in the number n of facets and the maximal number IfCv of source images for a facet. 

Proof. From the analysis in the proof of Theorem 16. 4[ the number of loop iterations 
is at most |src(f,S')| < |JF| srcy < n src^ . Within the main body of the algorithm, 
only standard geometric and linear algebraic operations are used, and these are all 
polynomial in n. Similarly, R,outine l6.2l uses only linear algebraic operations for every 
potential event E & S. Note that the cardinality of the set of potential events S during 
any iteration of the loop is bounded by |src(f , S*)! ■ < (^nsfc^) ■ = n^lfc^. 

Routine inini constructs Voronoi diagrams V(T) for finite sets T C M.'^. This com- 
putation requires |T|'^'^'^^ < ( srcy)'~'^'^^ time, which is polynomial for our fixed dimen- 
sion d. Therefore the total cost of the algorithm is also polynomial in n and WCy. □ 

Corollary 6.6. Let v and w be two points on the boundary S of the convex {d+ 1)- 
dimensional polyhedron P C R'^"'"^, and suppose that v lies interior to a facet. Then 
the geodesic distance ^{v,w) on S can be computed in time polynomial in n and IrCy. 

The restriction that v lie interior to a facet is unnecessary, and in fact Algorithm 16. II 
can be made to work for arbitrary points v\ see Sections 18.81 and 18.91 

Proof. Use Algorithm 16.11 to compute the foldout map ip : Uy ^ S. Find w' & % 
mapping to w = (p{w') G 5*, and compute the distance \v — w'\. By the isometry of 
the exponential map in Theorem 13.51 we conclude that ^{v,w) = \v — w'\. □ 

Remark 6.7. The complexity of Algorithm 16.11 is exponential in d if the dimen- 
sion is allowed to grow. For example, the number of vertices of P can be as large 
as n^^'^^ [Zie95j. Similarly, the number of cells in Voronoi diagrams of N points in MJ^ 
can be as large as N^^'^^ |Aur9H IFor95j . 
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On the other hand, for fixed dimension d Algorithm 16 . 1 1 can not be substantially im- 
proved, because the input and the output have costs bounded from below by (a poly- 
nomial in) n and src^, respectively. This is immediate for the input since P is defined 
by n hyperp lanes. For the output, we claim that the foldout in the output of Al- 
gorithm IHUl can not be presented at a smaller cost because it is a (usually nonconvex) 
polyhedron that has at least Ifc^ boundary ridges, meaning faces of dimension d — 1 
in the boundary of U^. To see why, let F be a facet with src^ source images, and 
for each i/ G srci? consider a shortest path whose sequential unfolding into Tp has 
endpoint v. If instead we sequentially unfold the paths into we get \stcf\ = sfCy 
segments emanating from v. Extend each of these segments to an infinite ray. Some 
of these infinite rays might pierce the boundary of through faces of dimension less 
than d—1, but adjusting their directions slightly ensures that each ray pierces the 
boundary of Uy through a boundary ridge. These ridges are all distinct because their 
corresponding rays traverse different facet sequences. 

Of course, the efficiency of Algorithm lti.ll does not necessarily imply that it yields an 
optimal solution to the discrete geodesic problem — or the unfolding problem, for that 
matter. (The problem of computing any nonoverlapping unfolding, not necessarily 
the source unfolding, is of independent interest in computational geometry [O'ROO^.) 
But although IfCy is not known to be polynomial in n, we conjecture in Section El 
that it is. See Section [8.101 for more history of the discrete geodesic problem. 

Remark 6.8. Following traditions in computational geometry, we have not specified 
our model of computation. In most computational geometry problems the model is 
actually irrelevant, since the algorithms are oblivious to it. In our case, however, the 
situation is more delicate, due to the fact that during each iteration of the loop we 
make a number of arithmetic operations that increase the error. More importantly, 
we make comparisons, which potentially require sharp precision. 

Let us note that Theorem 16 . 51 and its proof hold as stated for the complexity over M 
model [ BCSS98j . where there are no errors, and all arithmetic operations and com- 
parisons have unit cost. On the other hand, for the (usual) complexity over Z2 model 
jBCSSQSj, the details are less straightforward. Think of the hyperplanes as given 
by equations over the rational numbers Q, and suppose that the source point v is 
rational, as well. Then the vertices, ridges, source images, (tangents of angles in) 
angle sequences, and all other data throughout the algorithm are also rational. The 
whole computation can be done precisely over Q. The problem is keeping up with 
the computational precision as the denominators grow exponentially. Our preliminary 
calculations show that Theorem 16.51 still holds with several adjustments in the proof. 
Further exploration of this matter goes outside the scope of this paper and is left open. 

7. Convex polyhedral pseudomanifolds 

Recall the notion of polyhedral complex from Section ^ The results in Sections ^jHl 
hold with relatively little extra work for polyhedral complexes S that are substantially 
more general than boundaries of polytopes. Since the generality is desirable from the 
point of view of topology, we shall complete this extra work here. 
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Suppose that x is a point in a polyhedral complex S. Denote by 

Sxie) = {y e S \ fi{x,y) = e} 

the geodesic sphere in S at radius e from x. If (x) is the smallest face of S containing x, 
then for sufficiently small positive real numbers e, the intersection (x) fl Sx{£) is an 
honest (Euclidean) sphere {x)s of radius e around x. The set of points A^^; in 5* near x 
and equidistant from all points on {x)s is the normal space at x orthogonal to {x) in 
every face containing x. The spherical link of x at radius e is the set 

Nx{e) = {yeNx\f^{x,y) = €} 

of points in the normal space at distance e from x. When e is sufficiently small, the 
intersection of Nx{e) with any /c- dimensional face containing x is a sector inside a 
sphere of dimension k — 1 — dim(a;). The metric fi on S induces a subspace metric on 
the spherical link Nx{s). Always assume e is sufficiently small when Nx{e) is written. 

Definition 7.1. Let S" be a connected finite polyhedral cell complex of dimension d 
whose facets all have dimension d. Given a point x inside the union Sd-2 of all 
faces in S of dimension at most d — 2, we say that S is positively curved at x if the 
spherical link N^ie) is connected and has diameter less than ne. The space 5* is a 
convex^ polyhedral complex if S is positively curved at every point x G 5*^-2 ■ 

This is definition of positive curvature is derived from the one appearing in |Sto76j . 
It includes as special cases all boundaries of convex polyhedra; this is essentially the 
content of Proposition 11.21 

Spherical links give local information about geodesies, as noticed by D.A. Stone 
(but see also |BGP92[ Section 4.2.2]). 

Lemma 7.2 f |Sto76| Lemma 2.2]). Suppose S is a convex polyhedral complex. Then 
'J is a shortest path of length ae in the spherical link N^le) of a point x & S if and only 
if the union of all segments connecting points of 'j to x is isometric (with distances 
given by the metric on S) to a sector of angle a inside a disk in of radius e. 

Although Stone only uses simplicial complexes, we omit the straightforward gen- 
eralization to polyhedral complexes. Stone's lemma forces shortest paths to avoid 
low-dimensional faces in the presence of positive curvature. 

Proposition 7.3. Proposition^^ holds for convex polyhedral complexes S. 

Proof. Using notation from Lemma 17^ suppose that a < vr, and let 7 be the segment 
connecting the endpoints of 7 through the sector of angle a. Then 7 misses x. □ 

The rest of Section^ goes through without change for convex polyhedral complexes 
after we fix, once and for all, a tangent hyperplane Tp = M.'^ for each facet F. The 

■^Using 'convex' instead of 'positively curved' allows usage of the term 'nonconvex polyhedral 
complex' without ambiguity: 'nonpositively curved' is already established in the context of CAT(O) 
spaces to mean (for polyhedral manifolds, at least) that no point has positive sectional curvature in 
any direction. In contrast, 'nonconvex' means that some point has a negative sectional curvature. 
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choice of a tangent hyperplane is unique up to isometry. For convenience, we iden- 
tify F with an isometric copy in Tp, so that (for instance) we may speak as if F is 
contained inside Tp. This makes Definition II .51 in particular, work verbatim here. 

The main difficulty to overcome in the remainder of Sections ^jU] is the finiteness 
in Lemma 12.41 In the context of convex polyhedral complexes, this finiteness is 
fundamental. It comes down to the fact that shortest paths never wind arbitrarily 
many times around a single face inside of a fixed small neighborhood of a point. The 
statement of the upcoming Proposition 17.41 would be false if we allowed infinitely 
many facets, though it could still be made to hold in that case if the sizes of the 
facets and their dihedral angles were forced to be uniformly bounded away from zero. 

Proposition 7.4. Fix a real number r > and a convex polyhedral complex S. There 
is a fixed positive integer N = N{r, S) such that the facet sequence of each shortest 
path 7 of length r in S has size at most N . 

Proof. Pick a real number e > small enough so that the following holds. First, the 
sphere Sx{e) of radius e centered at each vertex x only intersects faces containing x. 
Then, for every point x on an edge but outside the union of the radius e balls around 
vertices, the sphere Sx{b/2) only intersects faces containing x. Iterating, for every 
point X on a face of dimension i but outside the union of all the previously constructed 
neighborhoods of smaller- dimensional faces, the sphere Sx{s/2'^) only intersects faces 
containing x. The existence of such a number e follows from the fact that every facet 
of S is convex, and that S has finitely many facets (Definition [7^. 

It suffices to prove the lemma with r = Let y be the midpoint of 7. The 

closed ball By{e/2'^~^^) of radius e/2'^~^^ centered at y intersects some collection of 
faces, and among these there is a face of minimal dimension k. Fix a point Xk lying in 
the intersection of this face with the ball By{e/2'^~^^). The ball Bxf.{s/2^) contains 7 
by the triangle inequality. However, Bx^{s/2^) might also contain a point Xj on a 
face of dimension j < k. If so, then choose j to be minimal. Iterating this procedure 
(at most d times) eventually results in a point x on face of dimension i such that 
Bx{e/2^) contains 7 and only intersects faces containing x. 

The metric geometry of 5* inside the ball Bx{e/2^) is the same as in Bx'{e/2^) for 
every point x' on the smallest face containing x, as long as Bx'{e/2^') only intersects 
faces containing x'. Since S has finitely many faces by Definition 17. !( we reduce to 
proving the lemma for shortest paths 7 after replacing S by the ball B = Bx{e/2^). 
In fact, we shall uniformly bound the number of facets traversed by any shortest 
path in B. For simplicity, infiate the metric by a constant factor so that B has radius 2. 
By a face of B we mean the intersection of B with a face of S. 

Note that B is isometric to a neighborhood of the apex on the boundary of a right 
circular cone when the dimension is d = 2. In this case, shortest paths in B can only 
pass at most once through each ray emanating from x. We conclude that the lemma 
holds in full (not just for B) when d = 2. Using induction on d, we shall assume that 
the lemma holds in full for convex polyhedral complexes of dimension at most d — 1. 

First suppose that x is not a vertex of S, so the smallest face {x) containing x 
has positive dimension. Then B is isometric to a neighborhood of x in the product 
(x) X A^^. of the face (x) with the normal space N^. Projecting 7 onto yields a 
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shortest path 7 whose facet sequence in the convex polyhedral complex A^^. has the 
same size as Induction on d completes the proof in this case. 

Now assume that a; is a vertex of S. If one of the endpoints of 7 is x itself, then 
7 is contained in some face of B. Hence we may assume from now on that x does 
not lie on 7. Consider the radial projection from B \ {x} to the unit sphere Sx{^) 
centered at x in B. If the image of 7 is a point, then again 7 lies in a single face; 
hence we may assume that radial projection induces a bijection from 7 to its image 
curve 7. Since the geometry of B is scale invariant, every path 7' in i? \ {x} mapping 
bijectively to 7 under radial projection has a well-defined facet sequence equal to C^. 

Choose another small real number e as in the first paragraph of the proof, but 
with B in place of S. Assume in addition that e < 1/2t[. Subdivide 7 into at least 
2*^1 e equal arcs, and use Lemma f7. 21 to connect the endpoints of each arc by straight 
segments in (the cone over 7 in) B. Lemma [7.21 implies that 7 has length at most vr, 
because 7 is a shortest path. Therefore each of the at least chords of 7 has 

length at most 2'^. The argument in the second paragraph of the proof now produces 
a new center x' for each chord, and we are assured that x' ^ x because the e-ball 
around x does not contain any of the chords. Hence the smallest face (x') containing x' 
has positive dimension, and we are done by induction on d as before. □ 

We shall see in Corollary 17.71 that Proposition 17.41 implies finiteness of the set of 
source images. But first, we need to introduce the class of polyhedral complexes for 
which the notion of source image — and hence the rest of Sections ^jHl — makes sense. 

Definition 7.5. A convex polyhedral complex S of dimension (i is a convex 'polyhedral 
pseudomanifold if S satisfies two additional pseudomanifold conditions: (i) each facet 
is a bounded polytope of dimension d, and (ii) each ridge lies in at most two facets. 

Remark 7.6. The "A.D. Aleksandrov spaces with curvature bounded below by 0" 
of [BGP92ij include convex polyhedral pseudomanifolds; see Example 2.9(6) there. 
Some of our results here, such as surjectivity of exponential maps and non-branching 
of geodesies, are general — and essentially local — properties of spaces with curvature 
bounded below by 0. But our focus is on decidedly global issues pertaining to the com- 
binatorial and polyhedral nature of convex polyhedral pseudomanifolds, rather than 
on a local analogy with Riemannian geometry. That being said, many of our results 
here can be extended to convex "polyhedral" pseudomanifolds with facets of constant 
positive curvature instead of curvature zero. We leave this extension to the reader. 

A flat point in an arbitrary convex polyhedral complex need not have a neighbor- 
hood isometric to an open subset of M*^, because more than two facets could meet 
there. In a convex polyhedral pseudomanifold, on the other hand, every fiat point not 
lying on the topological boundary has a neighborhood isometric to an open subset 
of M!^. This condition is necessary for even the most basic of our results to hold, 
including Corollary 12.21 (whose proof works verbatim for convex polyhedral pseu- 
domanifolds), and the definition of source image (which would require modification 
without it; see Section IHTSj) . 

We would have preferred to avoid the boundedness condition on facets, but the 
finiteness of the set of source images in Lemma IT^ can fail without it; see Section IH^ 
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Corollary 7.7. Lemma "KJ holds for convex polyhedral pseudomanifolds S. 

Proof. Since every facet is bounded, the lengths of all shortest paths in S are uniformly 
bounded. Proposition 17.41 therefore implies that there are only finitely many possible 
facet sequences among all shortest paths in S from the source. □ 

Corollary 17. 71 yields the following consequences, with the same proofs. 

Theorem 7.8. Proposition \2. 61 on the generalization of Mount's lemma and Theo- 
rem \2.yi on Voronoi diagrams hold verbatim for convex polyhedral pseudomanifolds S. 

The rest of Section |21 requires slight modification due to the fact that a convex 
polyhedral pseudomanifold S can have a nonempty topological boundary dS. 

Proposition 7.9. Fix a source point v in a convex polyhedral pseudomanifold S. 
Every warped point lies either in the topological boundary of S or in the cut locus . 

Proof. The same as Proposition 12. lOj assuming w is not in the boundary of S. □ 

In view of Proposition 17.91 the statement of Corollary 12. 1 II fails for convex polyhe- 
dral pseudomanifolds. Instead we get the following, with essentially the same proof. 

Corollary 7.10. Ifv is a source point in a convex polyhedral pseudomanifold S , then 

1. KyU dS is a polyhedral complex of dimension d — 1, and 

2. Ky U dS is the union Ky U Sd-2 U dS of the cut, warped, and boundary points. 

The considerations in Section El go through with one small modification: the non- 
compact flat Riemannian manifold 5*° is the complement in S of not just the {d — 2)- 
skeleton Sd-2, but also the topological boundary dS of S. The notion of what it 
means that a tangent vector at w & S can be exponentiated (Definition 13. 3|) remains 
unchanged, as long as w lies neither in 5*^-2 nor the boundary of 5*. Similarly, the 
notion of source interior ( Definition 13. 4|1 remains unchanged except that the exponen- 
tials exp(t() for < t < 1 must lie in neither the cut locus Ky nor the boundary dS. 

Theorem 7.11. Fix a source point v in the convex polyhedral pseudomanifold S . The 
exponential map exp : Uy ^ S on the source foldout is a polyhedral nonoverlapping 
foldout, and the boundary Uy^Uy maps onto KyUdS. Hence KyUdS is a cut set induc- 
ing a polyhedral nonoverlapping unfolding S^{Ky U dS) — > Uy to the source interior. 

Proof. Using Corollary 17. lUl in place of Corollary |2.1H the proof is the same as that of 
Theorem l3.5| except that every occurrence of S\Ky must be replaced by S\{KyUdS), 
and the open subspace 5"° must be defined as S* \ {Sd-2 U dS) instead of S" \ Sd-2- D 

Corollary 7.12. Every convex polyhedral pseudomanifold of dimension d is, as a 
metric space, obtained from a closed, star-shaped, polyhedral ball in by identifying 
pairs of isometric boundary components. 

Section m concerns local geometry in the context of convex polyhedra, and therefore 
requires no modification for convex pseudomanifolds, given that all of the earlier 
results in the paper hold in this more general context. 

In Section the only passage that does not seem to work verbatim for convex 
polyhedral pseudomanifolds is the proof of Proposition 15.51 That proof is presented 
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using language as if F and F' were embedded in the same Euclidean space R'^^^, as 
they are in the case S = dP. This embedding can be arranged in the general case 
here by choosing identifications of Tp and Tpi as subspaces of M'^"'"^ in such a way that 
the copies of F and F' intersect as they do in 5*. 

Finally, the algorithm in Section IHl works just as well for convex polyhedral pseu- 
domanifolds, as long as these spaces are presented in a manner that includes the 
structure of each facet as a polytope and the adjacency relations among facets. For 
example, folding maps along ridges shared by adjacent facets can be represented as 
linear transformations after assigning a vector space basis to each tangent hyperplane. 

For the record, let us summarize the previous three paragraphs. 

Theorem 7.13. The results in Sections^\^ and\^hold verbatim for convex polyhe- 
dral pseudomanifolds S in place of boundaries of convex polyhedra. 

8. Limitations, generalizations, and history 

The main results in this paper are more or less sharp, in the sense that further 
extension would make certain aspects of them false. In this section we make this 
sharpness precise, and also point out some alternative generalizations of our results 
that might hold with requisite modifications. Along the way, we provide more history. 

8.1. Polyhedral vs. Riemannian. The study of geodesies on convex surfaces, where 
d = 2, goes back to ancient times and has been revived by Newton and the Bernoulli 
brothers in modern times. The study of explicit constructions of geodesies on two 
dimensional polyhedral surfaces was initiated in |Lyu40| , and is perhaps much older. 

The idea of studying the exponential map on polyhedral surfaces goes back to 
Aleksandrov |Ale48t §9.5], who introduced it locally when d = 2. He referred to 
images of lines in the tangent space Tp to a facet F as quasi- geodesic lines on the 
surface, and proved some results on them specific to the dimension d = 2. Among his 
other results was the d = 2 case of Proposition 11.21 

A detailed analysis of the cut locus of 2-dimensional convex polyhedral surfaces was 
presented in |VP71j . This paper, seemingly overlooked in the West, gives a complete 
description of certain convex regions called 'peels' in |AU92j . which can be used to 
construct source unfoldings. The approach in |VP71j is inherently 2-dimensional and 
nonalgor it hmic . 

The study of exponential maps on Riemannian manifolds is classical jKob89j . 
F. Wolter |Wol85j proved properties of cut loci in the Riemannian context that are 
quite similar to our results describing the cut locus as the closure of the set of cut 
points. In fact, we could deduce part 2 of our Corollarv 12 . 1 II from |Wol85l Lemma 2] — 
in the manifold case, at least — using Proposition 11.21 fwhich has no analogue in Rie- 
mannian geometry). The method would be to "smooth out" the warped locus to 
make a sequence of complete Riemannian manifolds converging (as metric spaces) to 
the polyhedral complex S*, such that the complement of an ever decreasing neigh- 
borhood of the warped locus in S is isometric to the corresponding subset in the 
approximating manifold. Every shortest path to f in S" is eventually contained in the 
bulk complement of the smoothed neighborhood. 
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This method does not extend to the polyhedral case where S is allowed to be 
nonconvex, because Proposition 11.21 fails: shortest paths (between fiat points) can 
pass through warped points (Fig. E))- Moreover, the polyhedrality in the first part of 
Corollary 12 . 1 1 1 fails systematically when S is allowed to be nonconvex (see ^MMPST] ). 



Figure 14. Points in the shaded region have exactly two shortest 
paths to v; all of these paths go through the warped point w. 

8.2. Low-dimensional fiat faces. We assumed in Definition 17.11 that faces of di- 
mension d — 2 or less in convex polyhedral complexes must be nontrivially curved. 
Allowing convex polyhedral complexes where low- dimensional faces can be flat would 
break the notion of facet sequence in Corollarv ll.4| and would cause the set of warped 
points to differ from the union of all closed faces of dimension — 2, in general. The 
resulting definitions of folding map and sequential unfolding would be cumbersome if 
not completely opaque. Nonetheless, the resulting definitions would be possible, be- 
cause shortest paths would still enter facets (and in fact, all faces whose interiors are 
flat) at well-defined angles. The notion of exponential map would remain unchanged. 

Definition 12.31 and Theorem 12.91 should hold verbatim for the modified notion of 
convex polyhedral pseudomanifold in which low-dimensional fiat faces are allowed, 
because the generalized Mount Lemma (Proposition 12. 6p should remain true. Note 
that Mount's lemma relies mainly on Proposition 11.21 and Corollarv 12.21 The latter 
might be more difficult to verify in the presence of low- dimensional flat faces, because 
it needs every flat point to have a neighborhood isometric to an open subset of R*^. 
Thus one might have to assume S* is a manifold, and not just a pseudomanifold. 

Observe that Fig. IHl depends on not having low-dimensional flat faces: it uses the 
fact that the vertex bordering the shaded region must lie in the cut locus. 

8.3. Why the pseudomanifold conditions? Theorem 12.91 fails for convex poly- 
hedral complexes that are not pseudomanifolds, even when there are no flat faces of 
small dimension. Indeed, with the notion of cut point set forth in Deflnition l2.1l entire 
facets could consist of cut points. To see why, suppose there is a cut point interior 
to a ridge lying on the boundary of three or more facets, and note that the argument 
using Fig. 0] in the proof of Corollary 12.21 fails. For a more concrete construction in 
dimension d = 2, flnd a convex polyhedral pseudomanifold with a source point so that 
some edge in the cut locus connects two vertices (for example, take a unit cube with 




METRIC COMBINATORICS OF CONVEX POLYHEDRA 



37 



a source point in the center of a facet; see Fig. and then attach a triangle along 
that edge of the cut locus. The attached triangle ("dorsal fin") consists of cut points. 

The proof of Theorem 12 . 91 fails for non-pseudomanifolds 5* when we use the thinness 
of the cut set in the proof of Proposition The appropriate definition of cut point x 
for convex polyhedral complexes more general than pseudomanifolds should say that 
two shortest paths from x to the source leave x in different directions — that is, they 
pierce the geodesic sphere Sx{e) at different points. But Corollarv 12. 21 would still fail 
for shortest paths entering the "dorsal fin" constructed above. 

8.4. Aleksandrov unfoldings. The dimension d = 2 foldouts called 'star unfold- 
ings' in [ A092t .CH96, .AAOS97j were conceived of by Aleksandrov in |Ale48| p. 184]. 



Thus we propose here to use the term 'Aleksandrov unfolding' instead of 'star unfold- 
ing', since in any case these foldouts need not be star-shaped polygons. We remark 
that a footnote on the same page in |Ale48j indicates that Aleksandrov did not realize 
the nonoverlapping property, which was only established four decades later |AU92j . 

Aleksandrov unfoldings are defined for 3-dimensional polytopes P similarly to 
source unfoldings. The idea is again to fix a source point v, but then slice the 
boundary S of P open along each shortest path connecting f in 5* to a vertex. An 
example of the Aleksandrov unfolding of the cube is given in Fig. (see also Fig. EJ. 




Figure 15. Aleksandrov unfolding of the cube (the source point v is 
on the front face, while the left and back faces have no cuts). 

Note that when the source point is in the center of the face, the resulting Aleksandrov 
unfolding agrees with the source unfolding in Fig. ^ 

There is a formal connection between source and Aleksandrov unfoldings. Starting 
from the source unfolding, cut the star-shaped polygon f/„ into sectors — these are 
'peels' as in Section 18.11 — by slicing along the shortest paths to images of vertices. 
Rearranging the peels so that the various copies of v lie on the exterior cycle yields 
a nonoverlapping foldout |A()92j containing an isometric copy of the bulk of the cut 
locus. This rearrangement is illustrated in Fig. which continues Example 
(see |AAOS97] for further references). 

No obvious higher- dimensional analogue of the Aleksandrov unfolding exists, be- 
cause although the union of all shortest paths connecting the source point to warped 
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Figure 16. Source and Aleksandrov unfoldings of the cube, where the 
corresponding peels in both unfoldings are numbered from 1 to 8. 

points forms a polyhedral complex, this complex is not a cut set as per Definition 13 .11 
Indeed, thinking in terms of source foldouts again, the union of all rays passing from 
the origin through the images of warped points does not form the {d — l)-skeleton of 
a fan of polyhedral cones. Even when (i = 3, edges of S closer to the source point 
can make edges farther away look disconnected, as seen from the source point. An 
example of how this phenomenon looks from v is illustrated in Fig. El where the 
planes are meant to look tilted away from the viewer. 



Figure 17. Shortest paths to warped points making edges look disconnected. 

Circumventing the above failure of Aleksandrov unfoldings in high dimension would 
necessarily involve dealing with the fact that the set 3^,-2 of warped points generically 
intersects the cut locus in a polyhedral complex of dimension d — ?>. This "warped 
cut locus" usually contains points interior to maximal faces of making it impos- 
sible for these interiors of maximal cut faces to have neighborhoods in S isometric 
to open sets in M'^, even locally. Thus the picture in Fig. E[ where most of the cut 
locus can lie intact in is impossible in dimension (i > 3. The only remedy would 
be to make further slices across the interiors of the maximal faces of the cut locus 
before attempting to lay it flat in M"^. Making these extra slices in a canonical way, to 
generalize Aleksandrov unfoldings to arbitrary dimension, remains an open problem. 
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Figure 18. A source foldout with infinitely many source images 

8.5. Definition of source image. Some subtle geometry dictated our choice of 
definition of 'source image' fPefinition 12. 3^ . With no extra information to go on, 
we might alternatively have tried defining sicp as the (finite) set of endpoints of 
sequentially unfolded shortest paths 

• ending at a point interior to F; or 

• ending anywhere on F, including at a warped point. 

Both look reasonable enough; but the first fails to detect faces of dimension ci — 1 in 
the cut locus that lie entirely within ridges of S, while the second causes problems with 
verifying the generalized Mount Lemma (Proposition I2.(ij) as well as Proposition 15.51 
and Lemma 15.81 Not that the generalized Mount Lemma would be false with these 
"bonus" source images included, but the already delicate proof would fail. In addition, 
having these extra source images would add unnecessary bulk to the source poset. 

8.6. Finiteness of source images. As we saw in Lemma [2. 41 for boundaries of poly- 
hedra, or Proposition 17.41 and Corollary 17.71 for convex polyhedral pseudomanifolds, 
the number of source images is finite. The argument we gave in Lemma 12.41 relies 
on the embedding of S" as a polyhedral complex inside W^~^^ in such a way that each 
face is part of an affine subspace (i.e. not bent or folded). This embedding can be 
substituted by the more general condition that the polyhedral metric on each facet 
is induced by the metric on S (so pairs of points on a single facet are the same dis- 
tance apart in S as in the metric space consisting of the isolated facet). With this 
extra hypothesis, we would get finiteness of the set of source images even for convex 
polyhedral pseudomanifolds whose facets were allowed to be unbounded. However, 
allowing unbounded facets in arbitrary convex polyhedral pseudomanifolds can result 
in facets with infinitely many source images. 

For example, consider an infinite strip in the plane, subdivided into three substrips 
(one wide and two narrow, to make the picture clearer). Fix a distance i > 0, and glue 
each point on one (infinite) boundary edge of the strip to the point i units away from 
its closest neighhbor on the opposite (infinite) edge of the strip. What results is the 
cylinder 5* in Fig. ^| This cylinder would be a convex polyhedral manifold if its facets 
were bounded. The source foldout determined by a source point v in the middle 
of the wide substrip is depicted beneath S, shrunken vertically by a factor of about 2. 
The cut locus K^, which is a straight line along the spine of S, divides each substrip 
into infinitely many regions, so each substrip has infinitely many source images. 
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8.7. Generic source points. For generic choices of source point f , the source poset 
will be a chain — that is, a total order on events. The reason is that moving v infinites- 
imally changes differently the angles from different source images to ridges containing 
the same event point, precisely because these source images are sequentially unfolded 
along shortest paths leaving v in different directions. Note, however, that the dis- 
tances from these various source images to the same event point always remain equal. 

8.8. Warped source points. We assumed that the source point v & S lies in the 
relative interior of some facet; however, nothing really changes when v lies in the 
relative interior of some ridge. This can be seen by viewing the exponential map as 
living on the interior fiat points S° C S, as in Section IHl 

Moreover, simple modifications can generalize the exponential map to the case 
where v is warped. However, exponentiation on the complement of the cut locus 
cannot produce a nonoverlapping foldout in M'' if u e because the resulting cut 
locus would not be a cut set. Indeed, the cut locus would fail to contain all of 
so its complement could not possibly be isometric to an open subset of W^. On the 
other hand, exponentiation would instead produce a foldout of S onto the tangent 
cone to S at v. The main point is that the source point still connects to a dense set of 
points in S via shortest paths not passing through warped points, by Proposition II. 21 

8.9. Multiple source points. Let T = {vi. . . . ^Vk} C S* be a finite set of points 

on the boundary S = dP of a convex polyhedron P. Define the geodesic Voronoi 
diagram Vs{^) to be the subdivision of S whose closed cells are the sets 

Vsiy^Vi) = {w E S \ fi{vi,w) < nivj^w) for all 1 < j < k}. 

Just like the (usual) Voronoi diagrams, computing geodesic Voronoi diagrams is an 
important problem in computational geometry, with both theoretical and practical 
apphcations [AGSS89, KWR97. IFC98| (see |Mit()()l inn97j for additional references). 

Below we modify Algorithm 16.11 to compute the geodesic Voronoi diagrams in 5* 
when multiple source points are input. The modified algorithm outputs subdivisions 
of the facets of S that indicate which source point is closest. More importantly, it 
also computes which combinatorial type of geodesic gives a shortest path. In the code 
below, the while- DO loop and the routines remain completely unchanged. The only 
differences are in the initial and final stages of the pseudocode. 

Algorithm 8.1. (Computing the geodesic Voronoi diagram) 
INPUT convex polyhedron P C W^^^ of dimension + 1, and 
flat points f i, . . . , t>fc in the boundary S = dP 
OUTPUT geodesic Voronoi diagram Vs{T) in S 
INITIALIZE T^;' := {fj | Vi G F}, C^^^p '■= (F) for Vi e Tp, and 

£f '■= {ivi,F,R) G £f I POINT Vi e Tp can see R through F} 
[...] 

COMPUTE for each i: srcj := {{jy,F) | F G JF, z/ G Tp, and C,y^F begins with Fj}, and 
for each i: the subset Vs{T,Vi) := U{,.,F)esrc, ^C^f, z^) fl F of S* 
RETURN geodesic Voronoi diagram Vs{T) = (V5'(T, vi),. . . , Vs{T, Vk)) 
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That some of the source points Vi, . . . ,Vk might he in the same facet necessitates 
the caU to Routine 16.31 in the initiahzation of Sp- As we did before Theorem 16.51 
define Ifc to be the maximal number of source images for a single facet. 

Theorem 8.2. Let P C M.'^'^^ be a convex polyhedron and S = dP , with source points 
Vi, . . . ,Vk in S \ Sd-2- For fixed dimension d, Alaorithm \8.1\ computes the geodesic 
Voronoi diagram V5(T) in time polynomial in k, the number n of facets, and Ire. 

The proof is a straightforward extension of the proof of Theorems 16.41 and 16.51 it 
is omitted. Using observations in Section 18. 8t it is possible to modify Algorithm 18.11 
to work for set of arbitrary (that is, possibly warped) source points. 

8.10. The discrete geodesic problem. One of our motivating applications for this 
paper was to the discrete geodesic problem of computing geodesic distances and the 
shortest paths between points v and w in S. The reduction of this problem to com- 
puting source unfoldings is easy: construct the source foldout in the tangent cone 
at f , and compute the Euclidian distance between the images. 

We should mention here that for d = 2 essentially two methods are used in the lit- 
erature to resolve the discrete geodesic problem: the construction of nonoverlapping 
unfoldings as above (see |AAUS97| ICH96| ISS86j ). and the so-called 'continuous Di- 
jkstra' method, generalizing Dijkstra's classical al gorithm |Dij59| for finding shortest 
paths in graphs. The second method originated in |MMP87] and is applicable to non- 
convex surfaces (see also |Kap99| ). Interestingly, this method constructs an explicit 
geodesic wavefront, and then selects and performs 'events' one at a time. However, 
the time-ordering of events is based on the d = 2 fact that the wavefront intersects 
the union of ridges (edges, in this case) in a finite set of points. Our approach is a 
combination of these two algorithmic methods, which have previously been separated 
in the literature. We refer the reader to |MitOOj for more references and results on the 
complexity of discrete geodesic problems. In general, computing geodesic distances 
on arbitrary polyhedral complexes remains a challenging problem of both theoretical 
and practical interest. 

9. Open problems and complexity issues 

The source poset succeeds at time-ordering the events during wavefront expansion, 
but it fails to describe accurately how the wavefront bifurcates during expansion, 
because every event of radius less than r occurs before the first event of radius r 
in the source poset. On the other hand, the notion of 'geodesic precedence' from 
Definition 15.31 implies a combinatorial structure recording bifurcation exactly. 

Definition 9.1. Given a source point f on a convex polyhedral pseudomanifold S, 
the vistal tree T{v, S) is the set of events, partially ordered by geodesic precedence. 

The definition of geodesic precedence immediately implies that T(t>, S) is indeed a 
rooted tree. It records the facet adjacency graph of the polyhedral decomposition of 
the source foldout into cut cells of dimension d. Equivalently, this data describes 
the "vista" seen by an observer located at the source point — that is, how the visual 
field of the observer is locally subdivided by pieces of warped faces. Proposition 15.51 
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says precisely that the identity map on the set of events induces a poset map from 
the vistal tree to the source poset. In particular, when the source point is generic as 
in Section the source poset is a linear extension of the vistal tree. 

There are numerous interesting questions to ask about the vistal tree, owing to its 
geometric bearing on the nature of wavefront expansion on convex polyhedra. For 
example, its size, which is controlled by the extent of branching at each node, is 
important for reasons of computational complexity f Theorem 16. Sj) . 

Conjecture 9.2. The cardinality |src(f, S")! of the set of source images for a polyhe- 
dral boundary S is polynomial in the number of facets when the dimension d is fixed. 

Hence we conjecture that there is a fixed polynomial fd, independent of both S 
and u, such that |src(w, S) \ < fd{n) for all boundaries S = dP of convex polyhedra P 
of dimension d + 1 with n facets, and all source points v & S. Note that the cardinality 
in question is at most factorial in the number of facets: |src(f , S)\ < n ■ {n — 1)\ = n\. 
Indeed, each source image yields a facet sequence, and each of these has length at 
most n, starts at with facet F containing v, and does not repeat any facet. 

In fact, we believe a stronger statement than Coniecture 19.21 holds for boundaries 
of convex polyhedra. Given a shortest path 7, both of whose endpoints lie interior to 
facets, call the facet sequence traversed by 7 the combinatorial type of 7. 

Conjecture 9.3. The cardinality of the set of combinatorial types of shortest paths 
in the boundary S of a convex polyhedron is polynomial in the number of facets of S , 
when the dimension is fixed. 

That is, we do not require one endpoint to be fixed at the source point. The 
result is stronger than Conjecture 19.21 because source images are in bijection with 
combinatorial types of shortest paths in S with endpoint v. In both of the previous 
two conjectures, the degree of the polynomial will increase with d, probably linearly. 
When d = 2 both conjectures have been proved (see |AAQS97l ICH96j ) . 



The intuition for Coniecture 19.21 is that, as seen from the source point in a con- 
vex polyhedral boundary S, the faces of dimension d — 2 more or less subdivide the 
horizon into regions. (The horizon is simply the boundary of the source foldout f/„, 
as seen from v.) The phrase 'more or less' must be made precise, of course; and our 
inability to delete it altogether is a result of exactly the same phenomenon in Fig. El 
that breaks the notion of Aleksandrov unfoldings in higher dimension. 

The reason we believe Coniecture 19.31 is that we believe Coniecture 19.21 and there 
should not be too many combinatorial types of vistal trees. More precisely, moving 
the source point a little bit should not alter the combinatorics of the vistal tree, and 
there should not be more than polynomially many possible vistal trees. In fact, we 
believe a stronger, more geometric statement. It requires a new notion. 

Definition 9.4. Two source points are equivistal if their vistal trees are isomorphic, 
and corresponding nodes represent the same facet sequences. 

Again, the facet sequence corresponding to a node of the vistal tree is the list of 
facets traversed by any shortest path whose sequential unfolding yields the corre- 
sponding source image. Hence two source points are equivistal when their views of 
the horizon look combinatorially the same. 
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Conjecture 9.5. The equivalence relation induced by equivistality constitutes a con- 
vex polyhedral subdivision of the boundary S of any convex polyhedron. Moreover, the 
number of open regions in this subdivision is polynomial in the number of facets of S . 

Independent from the conjecture's validity, the vistal subdivision it speaks of — 
whether convex polyhedral or not — is completely canonical: it relies only on the 
metric structure of S. In addition, lower- dimensional strata of the vistal subdivision 
should reflect combinatorial transitions between neighboring isomorphism classes of 
vistal trees. Thus Conjecture 19.51 gets at the heart of a number of issues surrounding 
the interaction of the metric and combinatorial structures of convex polyhedra. 

Remark 9.6. An important motivation behind the above ideas lies in the computa- 
tion of the geodesic diameter of the boundary of a convex polytope. This is a classical 
problem in computational geometry, not unlike computing diameters of finite graphs 
(for the d = 2 case see |A092l lAAUSOTj ). One possibility, for example, would be to 
compute the vistal subdivision in Conjecture 19.51 and use this data to list the com- 
binatorial types of shortest paths. Each combinatorial type could then be checked to 
determine how long its corresponding shortest paths can be. Conjectures 19.31 and 19.51 
give hope that the geodesic diameter problem can be solved in polynomial time. 

Let us remark here that the polynomial complexity conjectures fail for nonconvex 
polyhedral manifolds of dimension d > 2. Note that this does not contradict the 
fact that when d = 2 there exists a polynomial time algorithm to solve the discrete 
geodesic problem (see section 15 . 101 above) . Indeed, the number of source images gives 
only a lower bound for our algorithm, while the problem is resolved by a different 
kind of algorithm. On the other hand, we show below that for d > 3 the discrete 
geodesic problem is NP-hard. The following result further underscores the difference 
between the convex and nonconvex case. 

Proposition 9.7. On (nonconvex) polyhedral manifolds, the number of combinatorial 
types of shortest paths can be exponential in the number of facets. In addition, the 
problem of finding a shortest path on a (nonconvex) polyhedral manifold is NP-hard. 

We present two proofs of the first part: one that is more explicit and works for 
all d > 2, and the other that is easy to modify to prove the second part. For the 
proof of the second part we construct a 3-dimensional polyhedral manifold, which is 
essentially due to Canny and Reif |CR87j . See Remark 19. 81 for comments on how to 
doctor these manifolds to make them compact and without boundary. 

Proof. To obtain a polyhedral domain with exponentially many shortest paths be- 
tween two points X and y, we consider a dimension d = 2 example. Simply take a 
pyramid shape polyhedral surface as shown in Fig. and observe that there exist 2'^ 
shortest paths between top point v and bottom vertex w, where k is the number of 
terraces in the pyramid. The omitted details are straightforward. 

Now consider a dimension d = 3 example of a different type. Polyhedrally sub- 
divide M.^ by taking the product of a line £ = M with the subdivision of in 
Fig. EHl Observe that there are only finitely many cells. Now add 4n hyperplanes 
Hq, . . . , H^n-i orthogonal to £, and equally spaced along C. This still leaves finitely 
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Figure 19. Nonconvex polyhedral surface in and shortest paths 
between points v and w. 




Figure 20. Polyhedral subdivision of a planar slice in 

many convex cells. Between hyperplanes H^k and -^4^+1, for all /c = . . . n — 1, remove 
all cells except the prisms whose bases are the top and bottom triangles in Fig. 1201 
Similarly, between hyperplanes and H^k+s, remove all cells except the prisms 

whose bases are the left and right triangles in Fig. 1201 

Now choose x and y to be points on i, with x being on one side of all the hyper- 
planes, and y being on the other side. Any shortest path connecting x to y must pass 
alternately through vertical and horizontal pairs of triangular prisms, and there is no 
preference for which of the two prisms in each pair the shortest path chooses. Thus 
the number of shortest paths is at least 4", while the number of cells is linear in n. 

For the second part, a construction in |CR87j presents a polyhedral domain B where 
the shortest path solution is NP-hard. This domain B is obtained by removing a set 
of parallel equilateral triangles from R^. To produce a manifold one has to thicken 
the squares into nearly flat triangular prisms. We omit the details. □ 

Remark 9.8. The polyhedral manifold S in the above proof is noncompact and has 
nonempty boundary; but with a little extra work, we could accomplish the same 
effect using a compact polyhedral manifold without boundary. The idea is to draw 
a large cube C around S in R'^, and place copies Ctop and Cbot of C as the top 
and bottom facets of a hollow hypercube inside R^. The remaining 6 facets of the 
hollow hypercube are to remain solid. The result is compact, but still has nonempty 
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boundary in Ctop and Cbot- This we fix by building tall 3- dimensional prisms in on 
the boundary faces, orthogonal to Ctop and Cbot, pointing away from the hypercube. 
Then we can cap off the prisms with copies of the cells originally excised from C C 
to get a nonconvex polyhedral 3-sphere in W^. 

Remark 9.9. The reader should not be surprised by the fact that computing the 
geodesic distance is NP-hard for nonconvex manifolds. On the contrary: in most 
situations the problem of computing the shortest distance is intractable, and in general 
is not in NP. For example, finding the shortest distance in a Cayley graph between 
two elements in a permutation group (presented by a list of generators in Sn) is 
known to be NP-hard even for abelian groups |EG81j . Furthermore, for directed 
Cayley graphs the problem is PSPACE-complete |Jer85j . 

Our final conjecture concerns the process of unfolding boundaries of convex polyhe- 
dra: if someone provides a polyhedral nonoverlapping foldout made of hinged wood, is 
it always possible to glue its corresponding edges together? Because wood is rigid, we 
need not only a nonoverlapping property on the foldout as it lies fiat on the ground, 
but also a nonintersecting property as we continuously fold it up to be glued. 

Viewing this process in reverse, can we continuously unfold the polyhedral bound- 
ary so that all dihedral angles monotonically increase, until the whole polyhedral 
boundary lies fiat on a hyperplane? This idea was inspired by recent works |BC02t 
ICDROSj and was suggested by Connelly.^ While the monotone increase of the dihe- 
dral angles may seem an unnecessary condition justified only by the aesthetics of the 
blooming, it is in fact crucial in the references above. 

As we have phrased things above, we asked for continuous unfolding of an arbitrary 
nonoverlapping foldout. But in fact, we only want to ask that there exist a foldout 
that can be continuously glued without self-intersection. Let us be more precise. 

Definition 9.10. Let S be the boundary of a convex polyhedron of dimension d+1 
in M'^+^. A continuous blooming of S" is a choice of nonoverlapping foldout U S, 
and a homotopy {0t : U — >• M.^~^^ | < t < 1} such that 

1. 00 is the foldout map U ^ S; 

2. (pi is the identity map on U ; 

3. 0t is an isometry from the interior f/ of f/ to its image, and 0t is linear on each 
component of the complement of the cut set in each facet, for < t < 1; and 

4. the dihedral angles between corresponding facets of (f)t{U) increase as t increases. 

An example of a continuous blooming is given in Figure |^ below. 

Conjecture 9.11. Every convex polyhedral boundary has a continuous blooming. 

Even though we ask only for existence, we believe that in fact the source unfolding 
can be continuously bloomed. As far as we know, this is open even for d = 2. 
Interestingly, we know of no nonoverlapping unfolding that cannot be continuously 
bloomed, and remain in disagreement on their potential existence. 



Private communication. 
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Figure 21. An example of a continuous blooming of the surface of the cube. 
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